THE CENTRAL LIMIT THEOREM FOR THE 
SMOLUCHOVSKI COAGULATION MODEL 



Abstract 



^ : Vassili N. Kolokoltsov* 

O 

<N 

bJ> 
< 

£N| , The general model of coagulation is considered. For basic classes of unbounded coagu- 

lation kernels the central limit theorem ( CLT) is obtained for the fluctuations around the 
p/j dynamic law of large numbers (LLN). A rather precise rate of convergence is given both 

PLh ■ for LLN and CLT. 

g 1 Introduction 

Throughout the paper we shall denote by X a locally compact topological space equipped 
j> ■ with its Borel sigma algebra and by E a given continuous non-negative function on X such 
0\ , that E(x) — > oo as x — > oo . Denoting by X° a one-point space and by X J the powers 
^|Ix...x X(j-times) considered with their product topologies, we shall denote by X their 
O \ disjoint union X = [JJ^qX^, which is again a locally compact space. In applications, X 
OO ■ specifies the state space of a single particle, X stands for the state space of a random number 
^ ■ of similar particles, and E describes some key parameter of a particle. In the standard model 
O ■ X — R + = {x > 0} and E(x) = x denotes the mass of a particle. 

J> ! By C(X) (respectively C^X)) we always denote the Banach space of continuous bounded 
[ functions on X (respectively its subspace of functions vanishing at infinity) with the sup-norm 
5_i 1 denoted by || • ||, by A4(X) - the Banach space of finite Borel measures on X with the norm also 
denoted by || • ||, and by M. + (X) - the set of its positive elements. The brackets (/, Y) denote 
the usual pairing (given by the integration) between functions / and measures Y, and \fi\ for 
a signed measure /i denotes its total variation measure. The elements of X will be denoted by 
bold letters, e.g. x = (xi, . . . , x n ) G X n C X. For a subset J in {1, . . . , n} we shall denote by 
|/| and J respectively its cardinality and its complement in {1, . . . ,n}, and by x/ the element 
of X' 7 ' given by the collection of Xi with i E I. 

Assume that we are given a continuous transition kernel K(x\,X2] dy) from X x X to X, i.e. 
a continuous function from X x X to M. + (X) (the latter equipped with its *-weak topology, 
i.e. the topology of the dual space to C 00 (X)). This kernel will be called the coagulation 
kernel and it will be assumed to preserve E, i.e. K(x\, x 2 ; dy) has support contained in the set 
{y : E(y) = E(x\) +E(x 2 )}. Moreover, K(xi,x 2 ; dy) is symmetric with respect to permutation 
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of xi and x 2 and has intensity K(xi,x 2 ) = j x K(xi, x 2 ; dy) enjoying the following additive 
upper bound: 

K{x 1 ,x 2 )<C{l + E{x 1 )+E{x 2 )) (1.1) 

with some constant C > and all X\,x 2 . 

The process of coagulation that we are going to analyse here is a Markov process Z(t) on 
X specified by the generator 

Lg(x) J2 I {9{xi,y)-9(x))K(xi;dy) (1.2) 

/C{l,...,n}:|J|=2 

(where x = (x±, . . . , x n )) of its Markov semigroup acting on an appropriate space of functions 
on X. It is known and not difficult to deduce from the theory of jump type processes (see 
e.g. [I]) that the process Z{t) is well defined by this generator (see e.g. a detailed probabilistic 
description of Z(t) in [28]). In the next Section the analytic properties of the Markov semigroup 
specified by L will be made precise. 
The transformation 

x = (xi, . . . ,x n ) h-> W x = h(S xl H S Xn ), (1.3) 

with h being a positive (scaling) parameter, maps X to the space A4hs{X) of positive measures 
on X of the form h5 x . By Zp we shall denote a Markov process on J\4hs(X) obtained from 
Z(t) by transformation (jl.3p combined with the scaling of L by h, i.e. Z^ is defined through 
the generator 

L h G g (h5^) = h Y / (^( x /' V) ~ 9(x))K{*i; d v) 

IC{l,...,n}--\I\=2 

= h Yl [(G g (h5* + h(5 y -8 Xl ))-G g (h5 x ))K(xr,dy) (1.4) 

i-C{l,...,n}:|/|=2 

on C(M hS (X)), where G g (h5 y ) = g(y) for any y G X. 

The law of large numbers dynamics (LLN) for the processes Z^ is given by the kinetic 
equation, whose most natural form is the weak one, i.e. it is the equation 

i(9^t) = l:[ [ (g(y) - g(x!) - g(x 2 ))K(x 1 ,x 2 ;dy)lJ>t( d Xi)lJ,t( d X2) (1.5) 

at 1 J XxX JX 

on fi t that has to hold for all g G C^X). It is known (see [28]) that if a family of initial 
measures hS^) f° r %t ls uniformly bounded with bounded moments of order (3 > 2, i.e. if 

sup / (1 + E^(y))hS x(h) (dy) < oo, (1.6) 

h Jx 

and if HS^qa tends *-weakly to a measure fio on X, as h — > 0, then the process Z^ with the 
initial data h6 x (h) tends weakly to a bounded solution /i t of (11.51) with initial condition /x that 
preserves and has bounded moments of order /3, i.e. such that 



sup 

s<t Jx 



[ {l + EP{y))».{dy)<oo (1.7) 
Jx 



and 

f E(y)fH(dy) = [ E(y)fi (dy) (1.8) 
J x J x 

for all t > 0. 

The first objective of this paper is to establish the corresponding central limit theorem 
(CLT), i.e. to show that the process 

Ft(Zj,to)=h- 1 '*{Z*{Z*)-K(ji )) 

of normalized fluctuations of Z^ around its dynamic law of large numbers fit converges in 
some sense to a generalized Gaussian Ornstein-Uhlenbeck process on M. (X) or a more general 
space of distributions. We obtain this result under some mild technical assumptions on the 
coagulation kernel thus presenting a solution to the problem 10 from the list of open problems 
on coagulation formulated in the well known review pQ. 

It is worth noting that though for the classical processes preserving the number of parti- 
cles (like interacting diffusions or Boltzmann type collisions) the results of CLT type are well 
established and widely presented in the literature (see e.g. jll] or [i5] and references therein), 
for the processes with a random number of particles the work on CLT began recently. For 
coagulation processes with discrete state space X = N and uniformly bounded intensities the 
central limit for fluctuations was obtained in [7j using stochastic calculus. For general processes 
of coagulation, fragmentation and collisions on X = R + , but again with bounded intensities, 
the central limit was proved by a different method in [19], namely by analytic methods of the 
theory of semigroups. The results of the present paper are obtained by developing further the 
approach from [19J. 

The second objective of the paper is to provide precise estimates of the error term both in 
LLN and CLT for a wide class of bounded and unbounded functionals on measures. Note that 
the usual "prove compactness in the Skorohod space and choose a converging subsequence" 
probabilistic method does not provide such estimates (see, however, [11] for a progress in this 
direction for interacting diffusions). Our main technical tool is the study of the derivatives of 
the solutions to kinetic equations with respect to initial data (this approach is inspired by the 
analysis of such derivatives for the Boltzmann equation in [IT]). The existence and regularity 
of these derivatives in weighted spaces of functions and measures are analyzed and the validity 
of CLT is proved to be connected with a certain kind of stability of these derivatives. The 
final estimates and their proofs depend on the structure and the regularity properties of the 
coagulation kernel. We demonstrate various aspects of our approach analyzing the following 
three classes of kernels: 

(CI) K(x x , x 2 ) = C(E(xi) + E(x 2 )). 

Remark. This is a warming up example, for the solutions to the main equations are given 
more or less explicitly in this case. 

(C2) K(x u x 2 ) < C(l + ^E\^)){1 + v^(^)). 

Remark. This model is analyzed to show the kind of results one can expect to obtain without 
assuming any differential or linear structure on the state space X. The unavoidable shortcoming 
of these results is connected with the absence of an appropriate space of generalized functions 
to work with. Hence the estimate of errors in LLN and CLT have to depend on something like 
the norm of F$ = (Zq — jj, )/Vh in Ai(X). But general /xo can not be approximated by Dirac 
measures Zq in such a way that Fq be bounded in A4(X). Hence the possibility to apply these 
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results beyond discrete supported initial measures /io is rather reduced. On the other hand, 
these kind of results are open to extensions to very general spaces. 

(C3) X = R + , K(xi,X2,dy) = K(xi, X2)8(y — x\ — X2), E(x) = x, K is non-decreasing in 
each argument 2-times continuously differentiable on (R+) 2 up to the boundary with all the 
first and second partial derivatives being bounded by a constant C. 

Remark. This is the case of our main interest. Unlike previous cases the estimate here 
turns out to depend on the norm of Fq coming from the dual space to continuously differen- 
tiable functions, and this norm can be easily made small for an arbitrary measure Hq on X. 
Therefore, to shorten the exposition, we shall prove CLT completely, up to the convergence 
of the distributions of processes on the Skorohod space of cadlag functions, only for this case, 
restricting the discussion of the first two cases only to the convergence of linear functionals. For 
simplicity, we choose here the state space X = R + of the standard Smoluchowski model, the 
extensions to finite-dimensional Euclidean spaces X being not difficult to obtain. Similarly we 
choose very strong assumptions on the derivatives (in particular, the kernels K(x, y) = x a + y a 
with a G (0, 1) are excluded by our assumption, as the derivatives of this K have a singularity 
at the origin). Finally let us stress that all kernels from (C1)-(C3) clearly satisfy fll.il) (possibly 
up to a constant multiplier). 

We refer to reviews [1] and [22J for a general background in coagulation models, and to [13] 
for simulation and numerical methods. 

The content of the paper is the following. In the next section we formulate the main results, 
and other sections are devoted to their proofs. In particular, Sections 4 and 5 are devoted 
to a detailed analysis of the equation in variations (linear approximation) around the solution 
of kinetic equation (JT3J) that describes the derivatives of the solution to (11.51) with respect to 
the initial measure fiQ. At the end of Sect. 5 a new property of the kinetic equation itself is 
established that is crucial to our proof of CLT, but seems to be also of independent interest. 
Namely Propositions 15. 5[ 15.81 show that the solution depends Lipschitz continuously on the 
initial measure in the topology of the dual to the weighted spaces of continuously differentiable 
functions or certain weighted Sobolev spaces. In Section 9 three general result are presented (on 
variational derivatives, on the linear transformation of Feller processes and on the dynamics of 
total variations of measures), used in our proofs and places separately in order not to interrupt 
the main line of arguments. In the Appendix some auxiliary facts on the evolutions specified 
by unbounded integral generators are presented. Though they should be essentially known to 
probabilists dealing with jump processes, the author did not find an appropriate reference. 

To conclude the introduction we shall fix the basic notations to be used throughout the 
paper without further reminder recalling as we go some relevant facts about Sobolev spaces 
and variational derivatives. 

(i) Weighted spaces of functions and measures arising from unbounded intensities of jumps. 
For a positive measurable function / on a topological space T we denote by Cf = C/(T) and 
Bf = Bf(T) (omitting T when no ambiguity may arise) the Banach spaces of continuous and 
measurable functions on T respectively having finite norm 

U\\f=U\\ Cf (T ) = snp(\<f ) (x)\/f(x)). 

x 

By C/ )00 = C/ )00 (T) and -B/,00 = Bf >00 (T) we denote the subspaces of C/ and Bf respectively 
consisting of functions such that (<p/f)(x) — ► as f(x) ^00. If / is a continuous function 
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on a locally compact space X such that f(x) — > oo, as x — > oo, then the dual space to C/ j00 (X) 
is given by the space AAf{X) of Radon measures on X with the norm \\Y\\f = sup{(0, Y) : 

W/<i}- 

We shall need also the weighted L p spaces. Namely, define L p j = L p j(T) as the space of 
measurable functions g on a measurable space T having finite norm ||<?||l p/ = ||<7//||l p - 

For X = R + = {x > 0} we shall use also smooth functions. For a positive / we denote 
by C) fi (X) the Banach space of continuously differentiable functions on X = R + such that 
lim x _> <f>(x) = and the norm 

IHIc£°(x) = WWcjix) 

is finite. By C 2 '°(X) we denote the space of two-times continuously differentiable functions 
such that lim x _> 4>( x ) — an d the norm 

H\\cy\x) = U'\\f + H"h 

is finite. By M.^{X) and Ai"j(X) we shall denote the Banach dual spaces to Cj'° and C 2 ' 
respectively. Actually we need only the topology they induce on (signed) measures so that for 
v E M(X)nM l f (X), i = 1,2, 

M\m).{X) =SU P {(0,Z/) : U\\ c yo (x) < 1}. 

Similarly one defines the spaces L*'° and L 2 '®, p > 1, as the spaces of absolutely continuous 
functions on X = R + such that lirn^o <t>{ x ) — with the norms respectively 

M\\l%w = U'Kax) = U'/f\\L p{ x), HWl^x) = H'/fh P (x) + WW/fi'Wux), 

as well as their dual (L 1 ^)' and (L 2 '°)'. 

As an example (needed later) let us estimate two of these norms for the Dirac measure 5 X 
on R + , x > and the function f(y) = f k (y) = 1 + y k : 

ll^ll^. fc (R + ) = sup{^ g(y) dy : \\g\\ Cfk < 1} = x + x k+1 /(k + 1); 



ll^ll(4;° fc )'(R + ) 



/X IPX 
g{y) dy : \\g/fk\\ L2 <i}<Jj Q fi(y) dy < c(k)^f k (x). (1.9) 

Not every v e M(X) belongs to Mj(X) or M 2 (X). Suppose that / is non-decreasing and 
v e M(X) is such that 



V\X) 



poo 

;)=/ i/(dj/) = o(l)(a;/(a;))- 1 , rr - oo. (1.10) 
Then by integration by parts for g G Cj' °(R+) 

/»oo /»oo 

(g,v) = - g{x)dv(x)= I g'{x)v(x)dx 
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(the boundary term vanish by (jl.lOp ). so that 

IIHI.Mi(x) 



// 



and 

Similarly, as 



\M){X) 



sup{0, v) 



<!}• 



/ 0(S)C/S< 

Jo 

it follows that if v G M.{X) is such that 



f q (y)dy 



i + i = i, 

v q 



5 = 0(1) 



-1/9 



oo, 



then 



1 1 
- + - 
V 1 



1. 



z/|| a 2,o y = sup{(^,^) : 



l p + II(^//) / I|l p <i}, p>i. 



In particular, recalling that the usual Sobolev Hilbert spaces H k (R) are defined as the comple- 
tion of the Schwarz space 5(R) with respect to the scalar product 



where 



(/, 9)m = (/, (1 - A) fc s) L2 = (1 + p 2 ) fe ^) L2 , 



denotes the usual Fourier transform, and that by duality (H k )' = H~ k it follows that 
IMI^oy = sup{(^,z>) : U/fWm < 1} = sup{(0,/z>) : \\<j>\\ m < 1} 



\\m 



H- 



\Hmip)\ 2 



dp 



1+p 2 



1.11) 



We shall use this formula in Section 7. 

(ii) Functional spaces describing indistinguishable particles. By C sym (X k ) we denote the 
Banach space of symmetric (with respect to all permutations of its arguments) continuous 
bounded functions on X k , and by C sym (X)- the Banach space of continuous bounded functions 
on X) whose restrictions on each X k belong to C sym (X k ). For a function / on X we denote 
by f® its natural lifting on X, i.e. f®(xx, . . . ,x n ) = f(xi) ■ ■ ■ f(x n ). 

If / is a positive function on X m = R™, we denote by Cj' sym (X m ) (respectively C^' sym (X m )) 
the space of symmetric continuous differentiable functions g on X m (respectively two-times 
continuously differentiable) vanishing whenever at least one argument vanishes, with the norm 



g || £yl,sym 



dg 



dxi 
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sup 



dg 



dXj 



(r 1 ) (x) 



and respectively 



\9\ 



c 



2.sym. 

f 



(X m ) 



dg 


+ 


d 2 g 


+ 


d 2 g 






C f {X m ) 


dx\ 


C f (X™) 


dx\dx2 


C f (X m ) 



(iii) Variational derivatives. For a function F on Aif(X) the variational derivative 5F is 
defined by 



5F(Y; x) 



lim 



-(F(Y + s5 x ) - F(Y)), 



where lim s ^o+ means the limit over positive s. Occasionally we shall omit the last argument 
here writing 5F(Y) instead of 8F(Y;.). The higher derivatives S F(Y; x%, x{) are defined 
inductively. 

As it follows from the definition, if 5F(Y; .) exists and depends continuously on Y in the 
*-weak topology of M. (or any JAf), then the function F(Y + s8 x ) of s G R+ has a continuous 
right derivative everywhere and hence is continuously different iable, which implies that 

F{Y + 5 X )-F(Y) = [ 5F(Y + s5 x ;x)ds. (1.12) 



We shall need an extension of this identity for more general measures in the place of the Dirac 
measure 5 X . To this end the following definitions turn out to be useful. For two continuous 
functions 4>, f such that < (p < f and f(x) — ► oo as x — > oo, we say that F belongs to 
C l {Mf t4) (X)), I = 0, 2, . . . , if F G C(M f ) and for all k = 1, I, S k F(Y; x x ,...,x k ) exists for all 
Xi, . . . ,Xk G X k , Y G A4f(X) and represents a continuous mapping M.f(X) \— > C^ m . 0( ^ 00 (-^ fc ), 
where A^/(X) is considered in its *-weak topology. We shall write shortly C l (A4f(X)) for 
C l (Ai fj(X)). All necessary formulae on the variational derivatives in these classes are collected 
in Lemma [9. 11 

Remark. The introduction of the cumbersome notations C m (A4f^(X)) is motivated by 
the fact that (under our assumption on the growth of the coagulation rates) if one considers 
the solution to the kinetic equations fi t with /x G A4 1+E p, then usually fi t G M.i + ep-^ and 
the derivatives of fit with respect to the initial data belong to C 1+ £fc with certain k < (3, see 
Sections 4 and 5. 

(iv) Propagators. If S t is a family of topological linear spaces, t G R + , we shall say that 
a family of continuous linear operators U t,r : S r i— > S 1 , r < t (respectively t < r) is a propa- 
gator (respectively a backward propagator), if [/*'* is the identity operator in S 1 for all t and 
the following propagator equation (called Chapman- Kolmogorov equation in the probabilistic 
context) holds for r < s <t (respectively for t < s < r): 

jjt^r = u t,v n q 13 ) 

By c and k we shall denote various constants indicating in brackets (when appropriate) the 
parameters on which they depend. 

For an operator U in a Banach space B we shall denote by ||Z7||b the norm of U as a bounded 
linear operator in B. 

At last, we shall use occasionally the obvious formula 



Jc{l r ..,n},[J[=2 ^ ^ 

valid for any / G C^ m (A 2 ) and x = (x 1; . . . , x n ) G A n . 



;i.i4) 
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2 Results 

First we recall some known results on the Cauchy problem for equation ( 11. 51) . A proof of the 
following two results can be found in [28J and p2] respectively. Recall that we always assume 
that our continuous coagulation kernel K{x\,X2] dy) preserves E and enjoys the estimate ( 11.11) . 

Proposition 2.1 If a finite measure /xq has a finite moment of order (3 > 2, i.e. if 

[ {1 + E^{y))fi {dy)<oo, (2.1) 
J x 

then equation ( 11.51) has a unique solution fi t with the initial condition /xo satisfying fl 1.7ft and 
(11.81) for arbitrary t. Moreover, 

sup/ E^{y)fi s {dy)<c{C,t,(3,{l + E,fi Q )){E^,fi Q ) (2.2) 

s<t J x 

with a constant c, and the mapping /xo fit is Lipschitz continuous so that 
sup||/x s (/xJ)-/x s (/x2)|| 1+BW < c{C,t,P,{l + E,fil + fil)){l + E 1+ ",fil + fil)\\fil-fi 2 \\ 1+E . (2.3) 



s<t 



for any to G [1, (3 — 1]. 

Proposition 2.2 Solutions fit from the previous Proposition enjoy the following regularity 
properties: 

(i) for any g G -Bi+^oo (respectively g G B 1+E p-i the function j g(x)fi t (dx) is a contin- 
uous function oft (respectively continuously differentiable function oft and (JT3J) holds); 

(ii) the function t i— > fi t is absolutely continuous in the norm topology of Ai 1+E p-i(X) and 
is continuously differentiable and satisfies the strong version of (11.51) in the norm topology of 

+Ef 3-j(X) for any 7 G {1,(3}. 

Remarks. 

1. The basic ideas of proving Proposition 12.11 go back to the analysis of the Boltzmann 
equation in [29]. Formulas (12.21) . (12. 3p are proved in [28] only for [3 — 2 and to — 1 
respectively, but the above extension is straightforward. 

2. Statement (ii) of Proposition 12 . 1 1 is proved in [T6] only for 7 = (3, but the extension given 
above is straightforward. In fact (ii) is done in the same way as the similar statement of 
Theorem IA.2I from Appendix. 

It is worth to observe that the operator has the form of the r.h.s. of equation (lA.lj) from 
the Appendix with Aihs instead of X and with the (time homogeneous) intensity 

a{h5x) = h f K{x I ;dy)<3Ch- 1 {l + E,h6 x ){l,h6 x ). (2.4) 

IC{l,...,n}:\I\=2 J 
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As the jumps in fll .4p increase neither (l,hS x )) nor (E,hS x ), it is convenient to consider the 
process Z\ on a reduced state space 

MIT 1 = {Y£ M hS : (1, Y) < e , (E, Y) < e x }. 

On this reduced space the intensity (12.41) is bounded (not uniformly in h). Hence L h is bounded 
in C(Ai e ^ ei ) and generates a strongly continuous semigroup of contractions there, which we 
shall denote by T/\ 

Let T t be a semigroup specified by the solution of ( 11.51) . i.e. T t /(/x) = f(^t), where \i t is the 
solution of fll.5p with the initial condition /i given by Proposition 12.11 with some (3 > 2. We 
can formulate now our first result. 

Theorem 2.1 [The rate of convergence in LLN] Let g be a continuous symmetric function on 
X m and F(Y) = (g,Y® m ). Assume Y = hS K belongs to -M^' ei , where x = (x\, . . . ,x n ). Then 
under the condition (CI) or (C2) 

sup\T t h F(Y)-T t F(Y)\ 

s<t 

< hK(C, m, k, t, e , ei)||(7|| ( i +B * )8 m(l + E k+3 , Y)(l + E 3 , Y)(l + E k , Y)" 1 " 1 (2.5) 
for any k > 1 and under the condition (C3) 
swp\T t h F(Y)-T t F(Y)\ 

s<t 

< hK{C, m, k, t, e , ei)||^|L2, SHm ™(1 + E k+ \ Y){1 + E k+ \ Y)(l + E k , Y) m ~ l (2.6) 

for any k > with a constant k. 
Remarks. 

1. We give the hierarchy of estimates for the error term making precise an intuitively clear 
fact that the power of growth of the polynomial functions on measures for which LLN 
can be established depends on the order of the finite moments of the initial measure. 

2. The estimates in case (C2) can be improved. However, not going into this detail allows 
one to keep unified formulae for cases (CI) and (C2). 

3. In Section 7 we prove the same estimates (12. 5|) . (12. 6p for more general functionals F (not 
necessarily polynomial). 

Recall that 

^(^, / i )=/ i - 1 / 2 (^(4)-^(/i )) 

is the process of the normalized fluctuations. The main goal of this paper is to prove that as 
h —>■ this process converges to the generalized Gaussian Ornstein-Uhlenbeck (OU) measure- 
valued process with the (non-homogeneous) generator 

k t F{Y) = l - I j j'(5F(Y),6 y -5 Zl -5 Z2 )K(z 1 ,z 2 ;dy)(Y(dz 1 )fi t (dz 2 )+fi t (dz 1 )Y(dz 2 )) 

4 r 



+ 7 / / / (5 2 F(Y), {S y - S Z1 - 5 Z2 r 2 )K(z u z 2 - dy) ^(dz,) ^ t (dz 2 ) . (2.7) 
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The generalized infinite dimensional Ornstein- Uhlenbeck processes and the corresponding 
Mehler semigroups represent a widely discussed topic in the current mathematical literature, 
see e.g. [21] and references therein for general theory, [27] for some properties of Gaussian 
Mehler semigroups and [B] for the connection with branching processes with immigration. The 
peculiarity of the process we are dealing with lies in its 'growing coefficients'. We shall analyze 
this process by the analytic tools developed in Sections 4 and 5. Let us start its discussion 
with an obvious observation that the polynomial functionals of the form F{Y) = (g,Y® m ), 
g G C sym (X m ), on measures are invariant under A t . In particular, for a linear functional 
F(Y) = (g,Y) ' 



A t F(Y) = -J J J (g(y)-g(z 1 )-g(z 2 ))K(z 1 ,z 2 ;dy)(Y(dz 1 )fi t (dz 2 )+ f i t (dz 1 )Y(dz 2 )). (2.8) 

Hence the evolution (in the inverse time) of the linear functionals specified by the equation 
F t = —A t F t , F t (Y) = (g t , Y) can be described by the equation 

g(z) = -A t g(z) = - J J (g(y) - g(x) - g(z))K(x, z; dy)fM t (dx) (2.9) 

on the coefficient functions gt (with some abuse of notation we denoted the action of A t on 
the coefficient functions again by A t ). Let U t,r be the backward propagator of this equation, 
i.e. the resolving operator of the Cauchy problem g = —A t g for t < r with a given g r . As 
we shall show in Propositions 15.21 - 15.44 the evolution U t,r is well defined in C 1+E k (X) in cases 
(C1)-(C2), and in C 2 '° Ek (X) in case (C3). 

Theorem 2.2 [CLT: convergence of linear functionals] Under condition (CI) or (C2) 
suv\n(g,F*(ZZ,fi )) - (U°' s g,F*)\ 



s<t 



< k(C, t, k, e , ex)Vh\\g\\ 1+E k(l + E k+ \ Z% + fi ) 2 ( 1 

for all k > 1, g G C 1+E k(X) , and under condition (C3) 
BuplE^.JffC^./io)) - (U°' s g,F h )\ 



Vh 



[2.10) 



M 1+ „ fc+1 (X), 



s<t 



<«(CUfc,e o ,e 1 )vft||0|| oa ,o ,(l + ^ fc+5 ,^ + ^o) 3 1 + 



1+E K 



Vh 



(2.11) 



M 1 



1+E k + 



for all k > 0, g G C 2 '® Ek (X), where the bald E denotes the expectation with respect to the 
process Z\ . 

To shorten the exposition, we shall deal in the future only with the most important case 
(C3). Though all the results have natural modifications in cases (CI) and (C2), let us stress 
again that for their applicability in cases (CI), (C2) one needs the initial fluctuation Fq to 
be bounded in the norm of A4i +E k+i(X), which is possible basically only for discrete initial 
distributions yUo- 
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For our purposes it will be enough to construct the propagator of the equation F = —A t F 
only on the set of cylinder functions C£ = C%{M™ +E k), tu = 1,2, on measures that have the 
form 

^ w "(y) = /((0i,y),...,(^,y)) (2.12) 

with / G C(R n ), and (pi, (p n G C™' k . By C fc we shall denote the union of for all n — 0,1, ... 
(of course, functions from C° are just constants). Similarly one defines the cylinder functions 
C%{{L™f +Ek )') under condition (C3). 

The Banach space of k times continuously differentiable functions on (with the norm 
being the maximum of the sup-norms of a function and all its partial derivative up to and 
including the order k) will be denoted, as usual, by C k (R d ). 

Theorem 2.3 [limiting Mehler propagator] Under the condition (C3) for any k > and a Ho 
such that (1+E k+1 , Ho) < oo there exists a propagator OU t,r of contractions onCk preserving the 
subspaces C k , n = 0,1, 2, ... such that OU t,r F, F G Ck, depends continuously on t in the topology 
of the uniform convergence on bounded subsets of M.™ +Ek , m = 1,2 (respectively (L™f +Ek )' in 

case k > 1/2) and solves the equation F = —A t F in the sense that if f G C 2 (R d ) in (12.121) . 
then 

4oC/*' r $t 1 ""^ n (y) = -AtOU*' r 0<t<r, (2.13) 



uni 



dt 

- m,0 



Iformly forY from bounded subsets of M.™ +Ek (respectively (L 2 ' x Ek )'). 

Our goal is to prove that this generalized infinite-dimensional Ornstein-Uhlenbeck (or Mehler) 
semigroup describes the limiting Gaussian distributions of the fluctuation process F/\ 

Theorem 2.4 [CLT: convergence of semigroups] Suppose k > and ho > are given such 
that 

sup(l + E fc+5 ,Z h + y u ) < oo. (2.14) 

h<ho 

(i) Let $ G C r k l (M\ +Ek ) be given by f T2TT2]) with f G C 3 (R n ) and all fa G C^ Ek {X). Then 
su V \E^(F t h (Z^,ii ))-OU ^(F^\ 



s<t 



< K(C,t, k,e ,ei)Vh max \\(pj\\ C 2,o ||/|| C 3 (R „)(1 + E k+5 , Z% + /i ) 3 1 



Vh 



(2.15) 



(ii) If $ G C^{M. 2 1+Ek ) (with not necessarily smooth f in the representation (I2.12p ) and F^ 
converges to some Fq as h — > in the -k-weak topology of M\ +Ek+1 , then 

lim \E^(F t h (Z^,Ho)) - Of/ '*<I>(F )| = (2.16) 

h— »0 

uniformly for Fq from a bounded subset of Ai\ +Ek+ i and t from a compact interval. 
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Theorem 2.5 [CLT: convergence of finite dimensional distributions] Suppose (12.141) holds, 
4>i,...,4> n G Cl'° Ek (R + ) and Ffr G {L%° +Ek+2 )' converges to some F in {L 2 ^ 1+Ek+2 )' , as h -> 0. 
Then the IL n -valued random variables 

®l...,t n = /i O )),... l (0„ l F*(Zf l /io))) l < t x < ... < t n , 

converge in distribution, as h —>■ 0, to a Gaussian random variable with the characteristic 
function 

n n pt n 

g tl ,...,t n (Pi,-,Pn) =exp{tY,Pj(U°' t ^ J ,F ) - Y, / ' £ PiPkH{s, U^fr, U s ^<j> k ) ds}, 

j=l j=l Jt i-l l,k=j 

(2.17) 

where to = and 



II(t, 0, V) = 7 / / /(0 ® V, («„ - - <^ 2 )® 2 )tf z 2 - dy) f i t (dz 1 ) l i t (dz 2 ). (2.18) 



4 

In particular, for t — t\ — ... = t n it implies 



n n n ft 

lim Eexp{i V(0„ F t h )} = exp{* V(£/°% F ) - V / ^'Vfc) 

j=l i=l j.^l" 70 

Theorem 2.6 [CLT: convergence of the process of fluctuations] Suppose the conditions of The 
' 2.J\ hold, (i) For any <p G C 2 ' fife (R + ) £/ie rea/ valued processes (0, F^(Zq, //q)) converge 



orem 



the sense of the distribution in the Skorohod space of cddldg functions ( equipped with its standard 
Ji-topology) to the Gaussian process with finite- dimensional distributions specified by Theorem 
Iff. 51 (nj The process of fluctuations F^{Zq,^q) converges in distributions on the Skorohod 
space of cddldg functions D([0, T}; (L 2 '® k+2 (R, + ))') (with Ji-topology), where (L 2 ^ £;fc+2 (R + )) / 
is considered in its weak topology, to a Gaussian process with finite- dimensional distributions 
specified by Theorem \2.5[ 



3 Calculations of generators 

From now on we denote by \i t = /x*(/xo) the solution to (jl.5p given by Proposition 12.11 with a 
f3 > 2. To begin with, let us extend the action of T/ 1 beyond the space C(M^ s ' ei ). 

Proposition 3.1 For any positive eo, e.\ and 1 < I < m the operator is bounded in the 
space Cr 1+E i ) \m(Jvl < j^ ei ) and defines a strongly continuous semigroup there (again denoted by 
T] 1 ) such that 

\\ T t\\c {x+El . )m {M e ^) <exp{c(C,m,/)eit}. (3.1) 

Proof. Let us show that 

L h F(Y)<c(C,m,l)e 1 F(Y) (3.2) 
for Y = hS x and F(Y) = (1 + E\ Y) m . One has 

L h F(Y) = h f[(l + E l ,Y + h(S y - 5 Xl )) m -(1 + E\ Y) m ]K{^- dy). 

IC{l,...,n}:\I\=2 J 
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As 

(1 + E\ h(5 y - 8 Xi - 8 X .)) < h[(E( Xi ) + E(xj)) 1 - E l ( Xl ) - E l ( Xj )] 

< hcWE'-^x^Eixj) + E(xi)E l -\xj)} 

and using the obvious inequality (a + b) m — a m < c(m)(a m ~ 1 b + b m ) one obtains 

L h F{Y) < hc(m, I) Yl [t 1 + E ^ yr^HE^ix^Eixj) + E(x l )E l ' 1 (x J )) 

/C{l,...,n}:|/|=2 

+ h m (E l - 1 (x i )E(x j ) + E(x i )E l ~ 1 (xj)) m ]K(^. I ; dy) 

<c(C,m,l) f J[(l + E l ,Y) m -\E l -\ Zl )E(z 2 ) + E{ Zl )E l -\z 2 )) 

+ h m - l (E l - 1 (z 1 )E(z 2 ) + E(z 1 )E l - 1 (z 2 )) m ](l + E( Zl ) + E(z 2 ))Y(dz 1 )Y(dz 2 ), 

where we used (11.14p . By symmetry it is enough to estimate the integral over the set where 
E(zi) > E(z 2 ). Consequently L h F(Y) does not exceed 

cj[(l + E\ Y) m - 1 E l - 1 (z l )E(z 2 ) + h m -\E l -\ Zl )E(z 2 )) m }(l + E(z l ))Y(dz 1 )Y(dz 2 ) 
< c(l + E\ Y) m (E, Y) + h m - x c J E m{l - 1)+1 { Zl )E m {z 2 )Y{d Zl )Y{dz 2 ). 

To prove (I3.2p it remains to show that the second term in the last expression can be estimated 
by its first term. This follows from the estimates: 

(E m , Y) = h s ^E m {x i ) < h {^E\xi)\ ml = h l ~ m/ \E l ,Y) m/ \ 
(E m( - l -V +1 ,Y) < h-\E mi - l - l \Y){E,Y) < h-" 1 ^- 1 ^ (E l ,Y) m ^~ 1 ^ (E,Y). 

Once (13.21) is proved it follows from (12 .4p that Lh is bounded in C^ +E i^m(J^A e ^ ei ), and (13.11) 
follows from Theorem IA.1I (or Proposition lA.lj) . 



The following statement is a straightforward extension of the previous one. 

Proposition 3.2 The statement of Proposition l37l\ remains true if instead of the space Cn + ^i w 
one takes a more general space C, 1+Eh \ mi r 1+E ij .yn. 3 ■■ 

Next we shall calculate the generator C of the deterministic semigroup T t and compare it 
with Lh- 

Proposition 3.3 (i) If F e C x {Mx +e p : i +e i3-i{X)), then 

jTtF^o) = jF{im) = £F(jh), (3.3) 

with 



2 



CF{Y) = \f [ (6F(Y;y)-6F(Y;x 1 )-6F(Y;x 2 ))K{x 1 ,x 2 ;dy)Y(dx 1 )Y(dx 2 ). (3.4) 

' X JXxX 
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(ii) If the variational derivative S 2 F(Y; x,y) exists for Y G M.^,^ and is a continuous 
function of three variables (Y taken in its -k-weak topology), then for any Y = /i<5 x 

L h F(Y) - CF(Y) = ~J J (SF(Y; y) - 26F(Y; z))K(z, z; dy)Y{dz) 

+ h 3 I (1 -s)ds V / (6 2 F(Y + sh(6 y - 6 Xl ); ■), (S y - 6^)® 2 )K(x i; dy). 

(3.5) 



/C{l,...,n}:|/|=2" 



(in) IfFe C(M{X)), Y = M x; then 

L h F(Y) = ±-J j J[F(Y + h(6 y -6 Zl -6 Z2 ))-F(Y)]K(z 1 ,z 2 ;dy)Y(dz 1 )Y(dz 2 ) 
1 

~~ 2 




[F(Y + h{5 y - 28 z )) - F{Y)}K{z, z; dy)Y(dz). (3.6) 



In particular, if F(Y) = (<fi,Y) with a continuous function (ft, then 

L h F(Y) = 1 [ [ [ [<f>{y) - </>(zi) - <P(z2)]K( Zl , z 2 ; dy)Y(dz l )Y(dz 2 ) 




[<j>(y)-2<j>(z)]K(z,z;dy)Y(dz). (3.7) 



Proof, (i) Follows from (19.31) and Proposition 12. 2( i). 
(ii) Applying (BT2|) (a) to yields 



LhF(Y) = h 2 Yl I A " <W^(x 7 ; dy) 

JC{l,...,n}:|J|=2 

+ h 3 f\l-s)ds j ^ F{ y + s/l ^ - ^ •' 0> - ^J® 2 )^; dy). 

^° /C{l,...,ni:|/|=2^ X 



Transforming the first term of the r.h.s. of this equation by (11.141) . yields (13. 5p . 
(hi) Is obtained by applying (11.141) directly to (jl.4p . 



Proposition 3.4 T/ie backward propagator 

fl 



U"r : C(fi P h (MS ,ei )) ^(AC/ 1 )) 



o/ t/ie process of fluctuations Fp obtained from by the deterministic linear transformation 
Q%(X) = h- l / 2 {Y - /2 t ), is given by 

U h f f> r F = (Q^T^n^F, (3.8) 

where VL^FiY) = F(Q!f:Y), and satisfies the equation 

dyhpsp = u^php. t<s<Tj (3.9) 
ds J 1 
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forF G C 3 (M i +E p^ +E i}-\{X)), where 

Vh 



A?F(y) = A t F(Y) + 



Vh 




(5F(Y), 5 y - 5 Z1 - 5 Z2 )K(z u z 2 - dy)Y(dz 1 )Y(dz 2 ) 




~\- I I (5F(Y),5 y -25 z )K(z,z;dy)(fi t + VhY)(dz) 

+ X ' ' I {6 2 F{Y),{5 y - 8 Z1 - S Z2 )^)K{z 1 ,z 2] dy){Y{dz 1 ) t x t {dz 2 ) + Y{dz 2 )^ 




h 
4 
h 
4 

Vh 





(5 2 F(Y), (6 y - 25 z f 2 )K(z, z; dy){m + VhY)(dz) 
(6 2 F(Y), (6 y - 5 Z1 - 6 Z2 f 2 )K( Zl , z 2 ; dy)Y (dzjY (dz 2 ) 



+ V±J(l-s) a dsJ J j(S 3 F(Y + sVh(5 y -5 Zl ~5 Z2 ),-),(5 y -5 Zl -5 Z2 r 3 ) 



'o 

x K(z 1 ,z 2 ;dy)(nt + VhY)(dz 1 )(fj, t + VhY)(dz 2 

h 3/2 



J (1 - s) 2 ds j j (6 3 F(Y + sVh(5 y - 24), •), (8 y - 26 z f 3 )K(z, z; dy){^ t + VhY){dz). 

(3.10) 



Proof. According to Lemma 19.21 the backward propagator Uff' r is given by (I3.8P and satisfies 
for F G C(fi[ 0j T](-A / l^' ei )) (see Lemma W7I\ for this notation), where 



Applying (13. 6p yields 



L h nMY) = Th 



Y + h(5 y - S Zl - 5 Z2 ) - /i t 



(3.11) 




Vh 



x K(zi, z 2 \ ; dy)Y(dz 1 )Y(dz 2 ) 




Y + h(5 y - 25 s ) - ^ 

Vh 



Y_-JH 
Vh 



Y-n t 
Vh 

K(z, z; dy)Y(dz) 



and consequently 



1 



M)- L L h Sl*F(Y) = - 




(F(Y+Vh(5 y -5 Z1 -5 Z2 )-F(Y) 
x K(zi, z 2 ; dy)(Vh~Y + fi t )(dz 1 )(VhY + /J, t )(dz 2 ) 



- l -J f[F(Y + Vh(S y -25 z )) - F(Y)]K(z,z;dy)(VhY + ii t )(dz). (3.12) 
Applying fl9~2D (b) yields 

F(Y + Vh(5 y - 5 Z1 - 5 Z2 )) - F(Y) = Vh(5F(Y), 5 y - 8 Zl - 5 Z2 ) + ^(5 2 F(F), (5 y - 5 Zl - 5 Z2 f 2 ) 



+ 



h 3/2 



(1 - s) 2 (5 3 F(Y + sVh{8 y - 5 Z1 - 5 Z2 )), (5 y - 5 Zl - 5 Z2 f 3 ) ds. 
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Hence developing the r.h.s. of f 1 3 . 1 2 j) in h yields the term at h x l 2 of the form 

\J J J (SF(Y),5 y - 5 Z1 - 5 Z2 )K(z 1 ,Z2]dy)fit(dz 1 )fi t (dz 2 ), 

the term at h° being precisely A t F(Y) given by (12. 7p . plus the remainder terms of order at 
least h}! 2 . As the above term of order h~ x l 2 cancels with the second term in (13.111) one obtains 

4 Derivatives with respect to initial data: existence 

This section is devoted to the analysis of the derivatives of the solutions to equation (11.51) with 
respect to the initial data. Namely we are going to study the signed measures defined as 

it = Ct(A*o! x; dz) = |^(yu ; x; dz) = lim -(^(/io + sS x ) - /it(yUo))- (4.1) 

0/iO s^0 + s 

We will occasionally omit some arguments in £ t to shorten the formulas. 

To motivate the formulation of rigorous results, let us start with a short formal calculations. 
Differentiating formally equation (II. 5p with respect to the initial measure Hq one obtains for £ t 
the equation 

-7r(fiS&)= / / {g{y) ~ g{x\) - g{x 2 ))K(x 1 ,X2]dy)i t {dx 1 )^ t {dx 2 ). (4.2) 
Jxxx Jx 

Of course, this is by no means a coincidence that this equation is dual to (12.91) . 
Introducing the second derivative 

r]t = T)t(%, w) = r] t (fi ; x, w; dz) = lim -(6(/i + sS w ; x) - &(/i ; x)), (4.3) 

s^0 + S 

and differentiating (14.21) formally one obtains for t] t the equation 
d 



,,(g,Vt(x,w; , ■)) = / / {g{y)-g{x 1 )-g(x 2 ))K(x x ,x 2 \dy) 

at J XxX JX 

x [rj t (x, w; dxi)fi t (dx 2 ) + £ t (x; dxi)i t {w] dx 2 )]. (4.4) 

The aim of this section is to justify these calculations and to obtain rough estimates for £ 4 and 
Vt- 

We start our analysis with a result on approximation of the solutions to kinetic equations 
by equations with bounded kernels. Let us introduce a cut-off kernel K n that enjoys the same 
properties as K and is such that K n [x\,x 2 \dy) = K(xi,x 2 ;dy) whenever E(xi) + E(x 2 ) < n 
and K n (x\, x 2 ) < Cn everywhere. 

For convenience, we shall assume (3 > 3 everywhere in this section. 

Proposition 4.1 Let /io i— > //" be the solution, given by Proposition ^. 1\ to the equation (11.51) 
with K n instead of K. Then /i™ — > fi t in the norm topology of A4i + e^(X) with w € [1,(3 — 1) 
and *-weakly in Ai 1+ EP uniformly for t from compact sets. 
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Proof. As the arguments given below use a rather standard trick in the theory of kinetic 
equations (similar ideas lead to a proof of Proposition 12. ip we shall give them only for u = 1. 

Let cr" denote the sign of the measure /i" — \i t (i.e. the equivalence class of the densities of 
/i" — \i t with respect to — ji t \ that equal ±1 respectively in positive and negative parts of 
the Hahn decomposition of this measure) so that — fi t \ = cr"(//" — /i 4 ). By Lemma one 
can choose a representative of a" (that we shall again denote by cr") in such a way that 

(1 + E, \fi - fi t \) = jf (<(1 + E), ±W - ^ ds. (4.5) 

Applying (11. 5p one obtains from (14. 5 j) that 

(1 + E, \£ -^\) = \f daj [«(1 + E)){y) - «(1 + E))^) - «(1 + E))(x 2 )] 

x [K n (xi, x 2 ; dy)fi™(dx 1 )n™(dx 2 ) - K(xi, x 2 \ dy)n s (dxi)fi s (dx 2 )]. (4.6) 

The expression in the last bracket in H4.6[) can be rewritten as 

(K n - K)(x 1 ,x 2 ; dy)^(dxi)(j n s (dx 2 ) 

+ K(x 1} x 2 ; dy)[((i"(dxi) - fj, s (dx 1 ))fj,^(dx 2 ) + ii s {dx 1 ){ii r ^{dx 2 ) - fi 8 (dx 2 ))]. (4.7) 

As \x n s are uniformly bounded in M. x+E ii and 

{1 + E{x 1 )+E{x 2 )) [ {K n -K){x 1 ,x 2 ;dy)<Cn- e {l + E{x 1 ) + E{x 2 )) 2+e 
Jx 

for 2 + e < (3, the contribution of the first term in (14.71) to the r.h.s. of H4.6[) tends to zero as 
n — > oo. The second and the third terms in (14.71) ar similar. Let us analyze the second term 
only. Its contribution to the r.h.s; of (14.61) can be written as 

\ f dsj [K(l + E)){y) - (of (1 + E)){ Xl ) - «(1 + E)){x 2 )\ 
x fT(xi,x 2 ; dy)o-™(xi)\fi n s (dxi) - ^ s (dxi)\^ n s (dx 2 ), 
which does not exceed 

\j d*j K 1 + E )iv) - (1 + *0(*i) + (! + ^Ofo)] 
x fT(a;i, x 2 ; dy)\n™(dx!) - (i a (dxi) \fj%(dx 2 ), 

because (cr™(xi)) 2 = 1 and |cr™(xj)| < 1, j = 1,2. Since K preserves E and (II .ip holds, the 
latter expression does not exceed 

C ds J(l + E(x 2 ))(l + E(xi) + E(x 2 ))\i^(dx 1 )- ^(dx^idx^ 

<C f dsil + E^-^Kh+E*. 
Jo 
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Consequently by Gronwall's lemma one concludes that 

||//" - fi t \\i+E = (1 + E, |/i" - /^ |) = 0(1)^00 exp <| t sup ||/i s ||i +jB 2 I . 

[ »6[0,t] J 

Finally, once the convergence in the norm topology of any J\A.x+e~i with 7 > is established, 
the *-weak convergence in Ai 1+E p follows from the uniform boundedness of \x n and /1 there. 

Proposition 4.2 (i) Under the assumptions of Proposition UTTl the backward propagator XJ t)T of 
equation ( 12.91) is well defined and is strongly continuous in the space Ci+s/s-i ooPO- Moreover, 
there exists a unique solution £ t to (14.21) in the sense that £0 = $x, & is a *-weakly continuous 
function {t > 0} 1— > M.x+ E i3-\(X) and ( 14.2ft holds for all g G Cx+e(X). Finally, 

\M",x)\\i+e» <K{t,\\vo\\i+E^){l + E"){x) (4.8) 

for all uj G [l,/3 — 1] and some constant k, £t is continuous with respect to t in the norm 
topology o/M 1+ e/3-i- f and is continuously differentiable in the norm topology of Jvlx+E0- 2 -^ for 
all e > 0. 

(ii) If £p are defined as £ f but from the cut-off kernels K n , then £™ — > £ t , as n — > 00 in the 
norm topology of M. l+E <^ with u G [l,/3 — 2) and in the *-weak topology of M. 1+Et3 -i. 

(Hi) £ t depends Lipschitz continuously on /i in the norm of Jvtx+E^ for uj G [1,(3 — 2] so 
that 

sup - M)\\i +E >* < ^(C,t,e , ei , (£ 2+w ,/4 + ^))||^ - f4\\i+^(l + E 1+ "(x)). 

8<t 

(iv) £ t can be defined by the r.h.s; of (14.11) with the limit existing in the norm topology of 
M i + E u (X) with uj G [1,(3 — 1) and in the *-weak topology of M-x+ep- 1 ■ 



Proof, (i) Equation (14.21) is dual to (12. 9p and is a particular case of equation (1A.14|) from 
Appendix with 



A t g(x) = / / (g(y) - g(x))K(z,x;dy)n t (dz), (4.9) 
Jx J x 



and 



B t g{x)= / g{z) I K(z,x;dy)nt(dz). (4.10) 
Jx Jx 

In the notations of Theorem IA.2I one has in our case 



a t (x) = / K(z,x;dy)n t (dz) < C(l + E(x))\\/j, t \\x+E, 
Jx Jx 



and for all u < (3 — 1 



< CIIsIIi+e- y (1 + B"W)(1 + )ft(dz) < 3C|| 9 || 1+E .|| ft || 1+£ . + .. 



Moreover, as uo > 1 



A t (l + E w )(x) <C f ((E(x) + E{z)f - E w {x)){l + E{z) + E(x))fi t (dz) 
Jx 

< Cc{uo) [ {E"~\x)E{z) + E?(z)){l + E(z) + E{x))^ t {dz) < Cc{u){\ + E w ){x)\\^ t \\ 1+E i + ^ 
Jx 

Hence the required well-posedness of the dual equations (12.91) and (14.21) and estimate (14.81) for 
uo =(3-1 follow from Theorem |Aj (i) , (ii) with ^i = 1 + E S , s £ [1,/? — 2), and ip 2 = l + E^ 1 . 
The last statement of (i) follows from Theorem I A. 21 (iii). Estimate (14.81) for other uo £ [1, (3 — 1] 
follows again from Theorem I A. 21 (i) and the estimates for at and B t given above. 

Remark. Note that 1 +E(x) should be of the order o(l) x ^ OQ (i' 2 1 ViX^) m order to fulfill the 
condition on the intensity a t from Theorem I A. 1[ Hence the necessity of the condition uo < (3 — 2. 

(ii) The proof is the same as the proof of Proposition 14.11 above. 

(iii) The proof of this statement is practically the same as for the corresponding statement 
(see Proposition I2.1( i)) for the solution of kinetic equation and uses the same trick as in the 
proof of Proposition 14. II above. Namely denoting £| = £t(/io),j = 1,2, one writes 



Ul - e t h +e - =J o ds J IMl + E«)){y) - (a s (l + E«))( Xl ) - (a s (l + EF))(x 2 )) 

x K(xi,x 2 ; dy) [£ s l s (dx l )n l s {dx 2 ) - ^(dx^n^dx^], 

where a s denotes the sign of the measure (again chosen according to Lemma l9T3il . Next, 

rewriting 

tl(dxi)vl(dx 2 ) - g(dxx)n 2 s (dx 2 ) = cr s (xi)\^ s - ^ 2 |(dxi)/i^(dx 2 ) + ^(dx^ifi] - £){dx 2 ) 

one estimates from above the contribution of the first term in the above expression for H^ 1 — 
£t 2 ||i+e- by 

j^ds j [E"(y) - E"{x x ) + E"(x 2 ) + l]K(x h x 2 ;dy)\e s ~ g^dx^idx*) 

< c(uo)C j ds j [Er-\ Xl )E{x 2 ) + E"{x 2 ) + 1](1 + E{x x ) + E(x 2 ))\g - ^\\(d Xl )^(dx 2 ) 

< «(C,o;,eo,ei) / ds\\£ - 

Jo 

and the contribution of the second term by 

n(C,uo,e ,ei) / ds\\f/ s - /j, 2 s \\ 1+ e^+i\\^\\i+e^ 
Jo 

< K,(C,uj,eo, ex, (E 2+UJ ,/j,l + f4))t\\nl - /i 2 J|i +i ^+i||£olli+i^+ 1 - 

It remains to apply Gronwall's lemma to complete the proof of statement (iii). 

(iv) General results on the derivatives of the evolution systems with respect to the initial 
data seem not to be applied directly for (II. 5ft . But they can be applied to the cut-off equations 



19 



(and this is the only reason for introducing these cut-offs in our exposition). Namely, as can 
be easily seen (this is a simplified "bounded coefficients" version of Proposition 2.2(h)), the 
solution /i" to the cut-off version of the kinetic equations (11. 5p satisfies this equation strongly 
in the norm topology of M. 1+E ri-^ for any e > 0. Moreover, ^ depends Lipshtiz continuously 
on hq in the same topology, the r.h.s. of the cut-off version of (11.51) is differentiate with respect 
to \i t in the same topology and £™ satisfies the equation in variation (14.21) in the same topology. 
Hence it follows from Proposition 6.5.3 of [23] that 

Q = £f G"o; x; dz) = lim -(//?(// , +sS x ) - M/^o)) 

s^0+ s 

in the norm topology of J\A 1+E/ 3- t with e > 0. Consequently 

{g, /4Vo + hS x )) - (g, Ht(no)) = / (g, 6?(^o + s5 x ; x; •)) ds 

Jo 

for all g G C 1+EI 3- eoo (X) and e > 0. Using statement (ii) and the dominated convergence 
theorem one deduces that 

(g,li t (fj, + h8 x )) - (g,/j, t (no)) = (g,£t(Ho + s6 x ;x; •)) ds (4.11) 

Jo 

for all g G Ci +E i(X) with 7 < (3 — 2. Again using the dominated convergence and the fact 
that £ t are bounded in M.\ +E ii-\ (as they are *-weak continuous there) one deduces that (14.111) 
holds for g G C 1+E /3-i iOC (X). Next, for these g the expression under the integral in the r.h.s. of 
(14. lip depends continuously on s due to Theorem IA.2I (iv) . which justifies the weak form of the 
limit (14. ip (in the *-weak topology of Aii +E /3-i). At last, by statement (hi) £t depends Lipshitz 
continuously on s in the r.h.s. of (14. lip in the norm topology of Aii +E ~t with 7 < (3 — 2. As £t 
are bounded in J\A 1+E /3~2 it implies that £ t depends continuously on s in the r.h.s. of (14.111) in 
the norm topology of A4i +E i with 7 < j3 — 1. Hence (14.111) implies (14. ip in the norm topology 
of M i +E i (X), 7 < f3 — 1, completing the proof of Proposition 14^21 

Proposition 4.3 (i) Under the assumptions of Proposition [Ol there exists a unique solution 
rjt to (14.41) in the sense that rjo = 0, f]t is a *-weakly continuous function t 1— > M.i +E (3-2 and 
(H~4D holds for g G C 1+E (X). Moreover 

\\rk(x,w; 0I|i+b» < n(C,t, \\fio\\ 1+ Ep) 

x sup {\\Ux; OHi+^+all^^l -)\\i+e + \\Uw; -)\\ 1+EU+a \\^(x] -)\\ 1+E ) (4.12) 

»e[o,t] 

/or 1 < c<j < /3 — 2 and some k. 

(ii) Ifrj™ are defined analogously to rjt but from the cut-off kernels K n , then 77™ — * rjt in the 
norm topology of M.i +E i with 7 < (3 — 3 and in the *-weak topology of j\4 1+E g-2 . 

(Hi) 7] t can be defined by the r.h.s. of (14 .3p in the norm topology of J\4\ +E i with 7 < (3 — 2 
and in the *-weak topology of A4 1+E/ 3-2. 

Proof, (i) Linear equation (14.41) differs from equation (14.21) by an additional non homoge- 
neous term. Hence one deduces from Proposition 12.11 (i) the well posedness of this equation 
and the explicit formula 

Vt {x,w)= / V t ' s Q s {x,w)ds, (4.13) 
Jo 
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where V t,s is a resolving operator to the Cauchy problem of equation (14. 2p given by Proposition 
I4.2( i) (or directly form Theorem IA.2j) and fl s (x,w) is the measure defined weakly as 



(g,Q s (x,w))= / (g(y) - g{x x ) - g(x 2 ))K(x 1 , x 2 ; dy)^ t (x; dx^^w; dx 2 ). (4.14) 

JXxX J X 

From this formula and the properties of £ t obtained above statement (i) follows. 

(ii) This follows from (14. 13j) and Proposition 12.2( h) . 

(iii) As in the proof of Proposition 12. 2( iv). we first prove the formula 

(g,£t(no + h5 w ;x, •)) - (g,£t(n ;x,-)) = / (g, r] t (fi + s5 w ; x, w; -))ds (4.15) 

Jo 

for g G Coo(X) by using the approximation 1]™, and the dominated convergence. Then the 
validity of (I4.15I) is extended to all g G Ci +E p-2 >00 using the dominated convergence and the 
above obtained bounds for rj t and £ t . By continuity of the expression under the integral in the 
r.h.s. of (I4.14p we justify the limit (I4.3P in the *-weak topology of J\4 1+E /3-2(X) completing the 
proof of Proposition I4.3I 



5 Derivatives with respect to initial data: estimates 

Straightforward application of Theorem lA.2l of the Appendix would give exponential dependence 
on (E",/j, ) of the constant k in (I4.8I) . And this is not sufficient for our purposes. The aim 
of this Section is to obtain more precise estimates for £ t . Unlike the rough results of the 
previous section that can be more or less straightforwardly extended to very general models 
with fragmentation, collision breakage and their non-binary versions (analyzed in [2], [15], [16]). 
the arguments of this section use more specific properties of the model under consideration. 

We shall use the notations of the previous section, assuming in particular that A t and B t 
are given by (14.91) . (I4.10p respectively. Due to the results of the previous section we are able 
to assume that all the Cauchy problems we are dealing with are well-posed. Recall that we 
denote by U t,r the backward propagator of the equation (12. 9p . 

Let us start with an estimate of the backward propagator U 1 / of the equation g = —A t g 
that holds without additional assumptions (C1)-(C3). 

Proposition 5.1 For all k > 0, U 1 / is a contraction in Cn +E k)-i and 

\U f /g(x)\ < K (C,fc,r,e ,e 1 )||^|| 1+i;fc [(l + J E fc )(a;) + ( J E fc+1 , / io)]. (5.1) 

Proof. U/ is a contraction in Cn +E k\-i by Proposition IA.ll because A t ((l + E k )~ l ) < 
(and this holds, because E k {y) > E k (x) in the support of the measure K(z,x; dy)). Next 

A t (l + E k )(x) < C J [(E(x) + E(z)) k - E k (x)}(l + E(x) + E(z))fi t (dz). 

Using the elementary inequality 

((a + bf - a k )(l + a + b)< c{k){a k {\ + b) + b k+1 + 1) 
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that is valid for all positive a, b, k with some constants c(k) yields 

A t (l + E k )(x) < Cc(k)[E k (x)(e + e x ) + e + (E k+ \ fi t )}. 

Then by (Q 

A(l + < k{C, k, t, e , e 1 )[E k (x) + 1 + (S fc+1 , /io)]- 

Hence (15.11) follows by Lemma IA.2I and the fact that U 1 / is a contraction. 

To simplify formulas we shall often use the following elementary inequalities : 

(a) (E\p)(E\v) <2(E k + l ~\v)(E,v), 

(b) (E k , v)E(x) < {E k+ \ v) + (E, v)E k (x) (5.2) 
valid for arbitrary positive v and k, I > 1. 

Proposition 5.2 Under condition (CI) suppose k>l. Then 

\U l ' r g{x)\ < K(C,k,r,e , ei )\\g\\ 1+E 4l + E k (x) + (E k+1 ,fi )(l + E(x))], (5.3) 
sup U t (jio;x, -)\\ 1+E * < k(C, t,eo,ei)[l + £ fe (x) + (1 + E(x))(E k+ \ /x )], (5.4) 

and 

sup||?7 s (/io;x,w;; < k(C, A;, t, e , ei) 

x [(1 + £ fe+1 (a;) + (^ +1 , /i )(l + ^ 2 (a;)) + (£ fe+3 , /i )(l + £(*)))(! + #H) 

+ (1 + £ fc+1 H + (E k+ \ + E 2 (w)) + (E k+ \ MQ )(1 + i?H))(l + E(x))\. (5.5) 

Proof. The simplicity of condition (CI) stems from the observation that the two dimensional 
functional space generated by the function E and constants is invariant under both A t and 
B t , and also the full image of B t belongs to this space. Hence representing the solution to 
g = -(A t - B t )g as 

9 = U^gr + 9 (5.6) 
one finds that g belongs to the above mentioned two dimensional space and satisfies the equation 

~g = -(A t -B t )g + B t U t /g r , g \ t=r = 0, t<r. (5.7) 

The corresponding homogeneous Cauchy problem 

= -(A t - B t )<f>, <j) r = a + (3E, 

can be written as 

a t + $tE(x) = Ca t (ei + (1, fi t )E(x)), a r = at,(3 r = (5 
in terms of (f) = a t + (3 t E(x) and clearly solves explicitly as 



P + a (l,fx s )e- ei{r - s) d& 



E(x), 
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which implies that 

||0t||i+.E < re(r,eo)||0 r || 1+ .E. 

It follows from (15. ip that 

\B t U\ r g r {x)\ <K(C,r,e ,e 1 )\\g r \\ 1+E 4B t (l+E k ) + (E k +\ii )B t l](x) 

< k{C, r, e , ei ) \\g r \\i+E*(l + {E k+ \ Mo ))(l + E{x)). (5.8) 

Solving the non-homogeneous equation (15. 7p by the Du Hamel principle and using the repre- 
sentation (I5.6P yields (15.31) . But by duality one gets 

Us(^o;x,-)\\ 1+E k =sup{(g,€ s (/i ;x,-)) : \\g\\i +E k < 1} 

= sup{(U°' s g,5 x ) : \\g\\ 1+E k < 1} = sup{U°' s g(x) : \\g\\ 1+E k < 1}, 

which implies (15. 4p . 
Now from ffl~T3l) 

sup||?7(/xo; x, w; -)||h-js* < ^ SU P w)\\ 1+E k =tsup sup (U s,t g, Q s (x, w)) 

s<t s<t s<t \g\<l+E k 

< k(C7, ^eo^Osupsupi^fi^tw)) : < 1 + E k (y) + (1 + E(y))(E k+1 , ^ )} 

< /<C,M,eo,ei)sup / f [1 + E h ( Xl ) + E h (x 2 ) + (E k+1 , fi )(l + E( Xl ) + E(x 2 ))] 

s<t J J 

(1 + E(xi) + E(x 2 ))£ s (x; dxi)£ s (w; dx 2 ). 

Dividing this integral into two parts with E{x\) > E(x 2 ) and E(xi) < E(x 2 ) one can estimate 
the first part as 

/tsup / / [l + E fc (x 1 ) + (E fc+1 ,/i )(l + J E(x 1 ))](l + J E;(x 1 ))e s (x;c/x 1 )es(w;rfx 2 ) 

s<< J J 

< Ksup -)ll (H&fo Olli+BH-i + (E fc+1 ,/i )||e s (x; .)||i +Jsa ) 



IKsl^ VII UlSs^i Vlli+Bfe+i t 



< k(1 + + £ fc+1 (x) + (1 + E(x))(E k+ \ /i ) 

+ (£ fe+ \ // )(1 + £ 2 (x) + (1 + £(x))(£ 3 , //„)] 

< «(1 + #H)[1 + E k+ \x) + (£ fc+1 , Mo )(l + E\x)) + (£ fe+3 , /i )(l + E(x))\, 

where we used both (15.21) (a) and (15.21) (b). As the integral over the second part is estimated 
similarly one arrives at (15.51) . 

Proposition 5.3 Under condition (C2) 

\\U^\\c 1+ ^ < exp{4C(t - r)(eo + ex)} (5.9) 

and the estimates (15.31) - (15.51) hold for all k > 1. 
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Proof. Since 

A t {l + y/E)(z) = J J (y/E(z) + E{x) - ^E{z))K{z, x; dy)fx t (dx) 

< C [ y/E(x)(l + v^M)(l + y/md)lh(dx) 
Jx 

< C(l + y^))(VE + E, fi t ) < C(e + 2ei)(l + y^E{zj), 

according to Proposition I A . 1 1 the positivity preserving backward propagator iFjf of the equation 
g = —A t g is bounded in C 1+v ^(X) with the norm not exceeding exp{C(t — r)(e + 2ei)}. On 
the other hand 

B t (l + Ve)(z) < C J (1 + v^M) 2 !! + \/E{z})fi t (dx) 

<2C(e + e 1 )(l + y/Ejz)). 

Hence B t are uniformly bonded in C 1+ ^(X) with the norm not exceeding 2C(eo + e\). Hence 
( 15.91) follows from the series representation (1A. 16|) for the backward propagator U r,t of the 
equation g = — (A t — B t )g. 

Now we use the same arguments as in the proof of Proposition 15.21 with 

\B t U\ r g r {x)\ < K (C,t- r,e^e 1 )\\g r \\ 1+& {l + y/E@j)(l + (E k+1 ,fi )) 

instead of (15.81) . Namely, in the representation of the solutions to g = —(A t — B t )g by the 
series ( 1A.16I) the first term is independent of B t and all other terms belong to C 1+ ^(X) and 
applying the above estimates for U*^ and B t in this space one deduces ( 15.31) . Other estimate 
follows now straightforwardly as in the previous Proposition (even with some improvements 
that we do not take into account). 

Proposition 5.4 Under condition (C3) for any k > the spaces C^ Ek and C 2 '® Ek (see Intro- 
duction for these notations ) are invariant under U t,r and 

(a) \(U t > r g)Xx)\< K (C,r,k,e ,e 1 )\\g\\ c i,o (1 + E k (x) + {E k+ \ //<,)), 

1+E k 

(b) \(U t > r g)"(x)\<K(C,r,k,e ,e 1 )\\g\\ G 2,o (1 + E k (x) + (E k+ \ /i )), (5.10) 

1+E k 

sup||£ s (/i ;x;-)IU </€(C,r,fc,eo,e 1 )[ J B(x)(l + (E fc+1 , / io)) + J B /£+1 (x)], (5.11) 

s<t 1+E 

and 

snp\\r] s (fx ;x,w; , -)\\ M 2 < K,(C,t, k } e , ei)(l + (E k+1 ,fi )) 

s<t ^+ Ek 

x {(E(x)(l + E k+2 , fio) + E k+2 (x))E(w) + (E(w)(l + E k+2 , /x ) + E k+2 (w))E(x)]. (5.12) 

Proof. Notice first that if g r {0) = 0, then g t = for all t according to the evolution 
described by the equation g = —(A t — B t )g. Hence the space of functions vanishing at the 
origin is invariant under this evolution. 
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Recall that E(x) = x in case (C3). Differentiating the equation g = —(At — B t )g with 
respect to the space variable x leads to the equation 

(g(x + z) - g(x) - g(z)) — (x, z)ii t (dz). (5.13) 

For functions g vanishing at the origin this can be rewritten as 

g'(x) = -A t g - D t g 

with 

DMx)= III rh(v)dv- I ch(v)dv) 

dx 



D t <fi(x)= I [ I (j)(y)dy- I 4>{y)dy ) ^(x,z)/i t (dz). 



o 



Since 

\\D t <p\\<2C\\<f>\\(E,v h ) = 2Ce 1 \ 
and U/ is a contraction, it follows from representation (1A.16[) with D t instead of B t that 



ll^' r |lcj'°(x) ^ K(C,r-t,eo, ei), 
proving (15.10[) (a) for k — 0. Next, for k > 

\D t (f>{x)\ < C\\<f>\\ 1+E * J {{x + z) k+l - x k+1 + z k+1 + 2z)fi t {dz) 
< Cc(k)\\(f>\\ 1+E u J (x k z + z k+1 + 2z)ii t (dz), 

which by (12.21) does not exceed 

c(C,A;,e 1 )||«/»|| 1+£ ,[(l+x fc ) + (£; fe+1 ) /.o)]. 

Hence by Proposition 15.11 

J* \U\ s D s U s /g{x)\ds < (r - t)K{C, r, k, e , ei) \\g\\ l+E u [1 + E k {x) + (E k+ \ //„)] , 
which by induction implies 

/ W^Ds, ■ ■ ■ D Sn U s /' r g(x)\d Sl ...ds n 

J t<si<---<s„<r 

( r — f) n 

<- T J -K n {C,r,k,e ,ei)\\g\\ 1+E k[l + E k (x) + {E k+ \^)\. 



Hence (15.10j) (a) follows from the representation (1A.16I) to the solution of (15. 13ft . 
Differentiating (15.131) leads to the equation 

g"(x) = -A t (g")(x)-iP t , (5.14) 
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where 

(g'(x + z) - g'(x)) — (x, z)ii t (dz) 

f ( f x+z f z \ d 2 K 

+ J \J 9 ^ dy ~J 9 (y)dy) -g^-(x,z)Vt(dz). 

We know already that for g r G C 2 +Ek the function g' belongs to 1 + E k with the bound given 
by (15. 101) (a). Hence by the Du Hamel principle the solution to (15.141) can be represented as 



As 

|Vt(z)l < K(C,r-t,e ,ei)(l + E k (x) + (E k+1 ,fi )), 

(15.1011 (b) follows, completing the proof of ( 15.101) . which by duality implies (15. lip . 
Next, arguing as in the proof of Proposition 15.21 one gets 

sup\\r] s (no;x,w;.)\\ M 2 < t sup sup{\(U s ^g, Q a (x, w))\ : \\g\\ c 2,o < 1} 

S<t i+E'* S<t 1+B fc 

< k(C, t, eo, ei) sup sup (g, Q s (x, w)), 
s<t gen fe 

where 

n fc = {g : g(0) = 0,max(|</(y)|, \g"{y)\) < 1 + E k (y) + (E k+ \^)}. 

It is convenient to introduce a two times continuously differentiable function \ on R such that 
x( x ) £ [0, 1] for all x, and x( x ) equals one or zero respectively for x > 1 and x < — 1. Then 
write Q s = Ql + Q 2 with Q 1 (respectively fl 2 ) being obtained by (14.141) with x( x i ~ x 2)K(xi,x 2 ) 
(respectively (1 — x( x i ~~ x 2))K(xi,x 2 )) instead of K(x±, x 2 ). If g G ITfc, one has 



((/,nj(a:,tw)) = / / [g{xx + x 2 ) - g(x l ) - g(x 2 ))x{xi~ x 2 )K(x x ,x 2 )^ s (x]dx l )^ s (w]dx 2 ), 




which is bounded in magnitude by 

d 

j - 

X2 



Us(w,-)\\ M i {x) SUp 



[(g(x! + x 2 ) - g{xi) - g(x 2 ))x(xi ~ x 2 )K(x 1 ,x 2 )]£ s ( x ; dxx) 



dx 2 

< Us{w,-)\\mI(x)Ms{ x ,-)\\m] , +1 (X) 



x sup 

Xl,X2 



d 2 



1 + E k+1 (xi)) 1 [(g(xi + x 2 ) - g(xi) - g(x 2 ))x(xi - x 2 )K( Xl , x 2 ) 

OX 2 OX\ 



Since 

d 2 



•[{g(x! + x 2 ) - g(x l ) - g(x 2 ))x{x x - x 2 )K(x 1 ,x 2 )] 
= g"{x x + x 2 )(xK)(xi,x 2 ) + (g'{xi + x 2 ) - g'(x 2 ))'- 



dx 2 dx\ 

d(xK)(xi,x 2 ) 



dx\ 

\ I 1 1 1 \ // ^d{xK){x u x 2 ) d 2 {xK)(x 1 ,x 2 ) 
+ (g (xi + x 2 )-g (an)) + (g(xi + x 2 ) - g(x x ) - g(x 2 ))- 



OX 2 OX\OX 2 
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this expression does not exceed in magnitude C(l + E k+1 (xi) + (E k+1 , /io)(l + -^(^i)) (up to a 
constant multiplier). Consequently 

1(^7,^(^,^)1 < ^(C7)||e £ (^, Oll^tiwII^C^' -)ll^ i+£7fc+1 (^)(l + (^^S/xo))- 

Of course, the norm of is estimated in the same way. Consequently (15. lip leads to (15.121) 
and completes the proof of Proposition 15.41 

We shall prove now the Lipschitz continuity of the solutions to our kinetic equation with 
respect to initial data in the norm-topology of the space M.\ +Ek - 

Proposition 5.5 Under the condition ( C3) for k > and m — 1, 2 

sup||/i a (^J) -fi 8 (j4)\\M m < ^(C } t,k } e 0} e 1 )(l + E 1+k ,^ + fiDW^l- ^ 2 \\ M m (5.15) 

s<t l+E 1+E 

Proof. By (JQ) and (JHUD 

{g,vM)~ tk(j%))= [ ds [ [ g(y)^ t (fi 2 + s(^ - fily,x;dy)(fil- n 2 )(dx). (5.16) 

Since 



o 



(<7,&(y;x;.)) = (£/ '*<7, &(*>;•)) = (t/ ^)(x), 

-im,0 



it follows from Proposition 15.41 that (g, £t(Y"; x; .)) belongs to C™? Fh as a function of x whenever 



(7 belongs to this space and that 

\\(g^ t (Y;x;.))\\ cm>0 (x < k(C, t, k, e , d) ||p|| cm ,o (x) (l + (£ fc+1 , Y)). 

Consequently (I5.15P follows from (15.161) . 

We shall discuss now the L 2 -version of our estimates. 

Proposition 5.6 Under condition ( C3) assume f is a positive either non-decreasing or bounded 
function. Then U £ are contractions in L 2j yf. (Thus U 1 / yield natural examples of sub- 
Markovian propagators with growing coefficients.) 

Proof. First observe that 

POO 

(u(x + y) — u{x)) g 2 (x)u(x)dx < (5-17) 



for any y > and a non-decreasing non-negative g (and any u, if only the integral is well 
defined). In fact, it is equivalent to 

(T y u, u) L21/g < {u,u) L ^ yg , 

where T y u(x) = u(x + y), which in turn follows (by Cauchy inequality) from (T y u, T y u) l 2 x , < 
(u,u)i 21/g . And the latter holds, because 

POO 

(T y u,T y u) L21/g = u 2 {x + y)g 2 (x)dx 
Jo 
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OO POO 

u 2 (z)g 2 (z — y) dz < / u 2 {x)g 2 {x) dx. 

Assume now that / is non-decreasing (and positive). From (I5.17P it follows that for arbitrary 
y>0 

poo 

(u(x + y) — u(x))K(x, y)f 2 (x)u(x) dx < (5.18) 



(we used here the assumed monotonicity of the kernel K), and hence (A t u, u)l 2 1/f < 0. Hence 
U 1 / can not increase the norm of L^i/f- To conclude that it is actually a semigroup of contrac- 
tions it remains to observe that due to Proposition 15.11 there exists a dense subspace in £2,1// 
that is invariant under U 1 ^ ■ Assume now that / is bounded. We again have to show the validity 
of (15.181) for a dense invariant subspace of functions u. First note that as the evolution U 1 / is 
well defined on continuous functions and preserves positivity and differentiability (by Proposi- 
tions 15.11 and 15 .4j) it is suffice to show (15.181) for positive functions u with bounded variation. 
Assuming that this is the case one can represent positive u as the difference u = u + — u~ of two 
positive non- decreasing functions (by decomposing its derivative in its positive and negative 
parts). As — (u~(x + y) — u~(x)) < 0, to show f 1 5 . 1 8 j) one needs to show that 



(u + (x + y) — u + (x))K(x, y) f 2 (x)u(x) dx < 0, 
and as — u~ is negative this in turns follows from 

POO 

I y = (u + (x + y) -u + (x))K(x,y)f(x)u + (x)dx < 0. 
Jo 

Denoting by M an upper bound for f 2 we can write 

POO POO 

ly = (u + (x+y)—u + (x))K(x,y)Mu(x)dx— (u + (x+y)—u + (x))K(x,y)(M—f 2 (x))u(x)dx, 



which is negative, because the integrand in the second term is positive and the first term is 
negative by (I5.17p . 

We shall consider now equation (5.13) that can be written in the form 

g'(x) = -A t (g')(x) - D\g' - D 2 g', (5.19) 

where 

PX / PX + Z fiTS \ 

(AV)(*) = y U <Ky)dy-^(x,z))K{dz), 

POO / PX + Z PZ \ QJ£ 

(D 2 (p)(x) = J \l x<z J (f)(y) dy- J <f)(y) dy\ ~^r(x, z)fx t (dz). 

Proposition 5.7 Under condition (C3) for any fk(x) = 1 + x k , k > 1/2, the spaces L™f k , 
ra— 1,2 (see introduction for these notations), are invariant under U t,r and 

WU^W^o </t(C,r,A;,e ,e 1 )(l + ( J E; fc+1 / 2 ,/io)), m = l,2. 

2 'fk 
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Moreover, for g G L™'® one can represent (t/* ,r g)' as the sum of a function from L™'® k with the 
norm not exceeding n(C, r, k, eo, ei)\\g\\ L m,o and a uniformly bounded function with the sup-norm 

not exceeding k(C, r, k, eo, ei)||g|| L m,o (1 + (E k+1 / 2 , fio))- Consequently 

2 'fk 

8up\\Uiio;x;-)\\ (L i,o y < K (C,r,fc,e ,e 1 )[E(x)( J E; fc+1 , / i ) + (l + J E; fe+1 / 2 (x)]. (5.20) 

s<t V 2 'V 

Proof. Let us show first that 

\\Dl\\ L2Jk <ex2 fe a (5.21) 
In fact, for a continuous positive <fi and an arbitrary z > 

pX~\~Z nit . * 2 

l|i*<* / <P(y) dy\\l 2 f = lim / uV^ + -)^ + -)Vt" 2 W (il 

Jx J i,j=l 

/n . . 

l z<x T(<P/fk)(x + 3 -)(<P/fk)(x + -)-2 2k dx, 
~ n n n z 

because 

1 2 k 2 k 

< < 

fk(x) fk{2x) f k (x+jz/n) 

for all j < n, z < x. Taking now into account that 

/n 
i 2 <, &/fk)(x + -)(Hfk)(x + -)dx< u/f k \\l 2 

one deduces that 

fX + Z 



\U<x I mdy\\l 2Jk <z^ k U\\l 2Jk . 



X 



Consequently 

IIAVIUv, < C / ||1 2 <, / 4>{y) dy\\ LW i t {dz) < C2 k e 

JO Jx 



2,/fc' 



which implies (I5.2ip . 

As clearly the same bounds hold for \\Dl ||c(r+)' the equation 



g'(x) = -A t (g')(x)-D}g' 

specifies a propagator C/* ,r , t G [0,r], of bounded operators in both C(R + ) and L 2 ,/ fc (R+) with 
uniform bounds depending on r, k, eo, e%. Next 

\{D 2 t <P){x)\ <2C j £ " \<P{y)\dyUdz) 
for all x, which by Cauchy-Schwartz inequality does not exceed 



2C|Mk* / \l r fl{y)dy^{dz)<CcmHL 2 Al + {E k+l l\^)). 
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Hence writing the solution to the Cauchy problem for equation 15.131 with <p r e L 2 j k as a 
perturbation series (1A.16|) with respect to perturbation D 2 one represents the solution as the 



sum <p\ + 02 with 

||0ilU 2 , /fe < c(/c,e o ,ei,r)||0 r || L2Jfc 

and 

||4l|c(R + ) < (l + (i? fc+1/2 ,/io))^(C,A;,e o ,e 1 ,r)||0 r || i2i/fc 

so that 

Hi\\L 2Jk < (l + (E k+1 / 2 ^ ))K(C,k,e ,e u r)U r \\ L2fk 
whenever k > 1/2. In particular for these k 

U^'l^o < (l + (E fc+1 / 2 , / i )) K (C,fc,eo,e 1 ,r). 

As (^ s (/io; ^; •)> flO = (^a;, ^°' s (7), this implies (I5.20p by (II. 9p . The evolution U t,r in the space 
L 2 2°f k is analyzed quite similarly. 

We conclude with the following analog of Proposition I5.5[ whose proof follows from Propo- 
sition 15.71 by the same argument as Proposition 15.51 follows from Proposition 15.41 

Proposition 5.8 Under the condition (C3) for k > 1/2 and m = 1,2 

sup ||/x s Oo) - Vs(l4) || (L -.o v < k(C, t, k, e , ei) (1 + E 1+k/2 , fi] + /i 2 ) - l4\\ (L ™>° y (5-22) 

6 The rate of convergence in the LLN 

Proof of Theorem \2.1\ Recall that ^t{Y) means the solution to equation (11. 5p with initial data 
Hq = Y given by Proposition 12.11 with a (3 > 2. We shall write shortly Y t = ^t{Y) so that 
T t F(Y) = F(jh(Y)) = FiY t ). For a function F(Y) = (g,Y® m ) with g e C ( ^ )0oo (X m ), 
m > 1, and Y = h5 x one has 

ft 

-\h T7T (\r\ I rph 



T t F(Y)-T t h F(Y)= / T t h L s {L h -C)T s F{Y)ds. (6.1) 
As T t F(Y) = (g,Y t ® m ), Propositions SJ and yield 

ST t F(Y;x) = m \ g(y h y 2 , . . . ,y Tn )£ t (Y;x;dy 1 )Y t ® (m ~ 1 \dy 2 - ■ -dy m ), 



and 

5(m— 1) 



5 T t F(Y;x,w) = ml g(yi,y 2 , ■ ■ ■ ,y m )vt(Y; x,w; dyi)Y t (dy 2 ■ ■ ■ dy m ) 
J x™ 

+ m(m-l) / g(y 1 ,y 2 ,...,y m )£ t (Y]X]dy 1 )£ t (Y]w-,dy 2 )Y® (m ~ 2 \dy 3 ---dy m ). (6.2) 



X m 



Let us estimate the difference [Lh — £)T t F(Y) using ( 13.51) (with T t F instead of F). Let us 
analyze only the more weird second term in ( 13.51) . as the first one is analyzed similar, but much 
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simpler. We are going to estimate separately the contribution to the last term of (13. 5p of the 
first and second term in (16. 2p . 

Assume that the condition (CI) or (C2) holds and a k > 1 is chosen. Note that the norm 
and the first moment (E, •) of Y + sh(S y — 5 Xi —8 X .) do not exceed respectively the norm and the 
first moment of Y. Moreover, for s G [0, 1], h > and Xj,y 6 X with E(y) = E(xj) + E(xj) 
one has 

{E\ Y + - 5 Xi - 6 X .)) = (E k , Y) + sh(E( Xi ) + E( Xj )) k - hE k (x t ) - hE k ( Xj ) 
< {E k ,Y) + hc(k)(E k - l (x i )E(x j ) +E{x i )E k - 1 {x j )) 

with a constant c(k) depending only on k. Consequently by Proposition 15.21 

\\r] t (Y + sh(5 y - 5 Xi - 5 Xj );x,w; < «(C, A;, t, e , e x ) (1 + E(w)) 

{1 + E fc+1 (x) + [(E k+1 , Y) + hc{k){E k {xi)E{xj) + E(x i )E fc (x j ))](l + ^ 2 (x)) 

+ [(E fc+3 , Y) + hcik^E^ix^Eixj) + Eix^E^ix^l + E{x))} + 

where by dots is denoted the similar term with x and w interchanging their places. Hence the 
contribution to the last term of (13. 5p of the first term in (16. 2p does not exceed 

k{C, t, k, m, e , ei)||0|| (1+E *)Bm(l + E k , Y^h 3 ^(1 + E( Xi ) + E( Xj )) 2 {l + E k+l { Xi ) + E k+1 { Xj ) 

+ [(E k+ \ Y) + hc(k)(E k (x t )E( Xj ) + Eix^E'ixjMl + E 2 ( Xi ) + E 2 {x 3 )) 
+ [(E k+3 , Y) + hcik^E^ix^Eixj) + E(xi)E k+2 (xj)](l + E( Xi ) + E( Xj ))}. 

Dividing this sum into two parts, where E(x{) > E(xj) and respectively vice versa, and noting 
that by the symmetry it is enough to estimate only the first part, allows to estimate the 
contribution to the last term of (13. 5p of the first term from (16.21) by 

k{C, t, k, m, e , ei) || g\\ {1+E k^m(l + E k , Y) m ~ l h 3 

{1 + E k+3 (x t ) + (1 + E\ Xi ))[{E k+ \ Y) + hc{k)E k {x t )E{x 3 )} 
+ (1 + E 3 (xi))[(E k+3 , Y) + hc(k)E k+2 (xi)E(xj)]}. 

The main term in this expression (obtained by ignoring the terms with hc(k)) is estimated by 

4g\\(i+E*)®m(l + E k , ^-"^[(l + E k+3 , Y) + (1 + E\ Y){E k+ \ Y) + (1 + E 3 , Y){E k+3 , Y)], 

where the first two terms in the square bracket can be estimated by the last one, because 

(E\Y)(E k+ \Y) < 2(E 2 ,Y)(E k+3 ,Y). 

It remains to observe that the terms with hc{k) are actually subject to the same bound, as for 
instance 

h^^i^Eix^il+E^Xi)) < h 2 (E k + E k+ \Y)(E,Y) < c(k)h(E k+3 ,Y)(E,Y) 2 . 
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Consequently, the the contribution to the last term of (13. 5p of the first term in (16. 2\i does not 
exceed 

hn{C, t, k, m, e , ei)||</|| (1+£ . fc)9m (l + E k , K) m_1 (l + E k+3 , Y){1 + E 3 , Y). (6.3) 
Turning to the contribution of the second term from (16. 2p observe that again by Proposition 



||f t (y + sh(5 y - 5 Xi - S x .);x;-)\\ 1+E k < k(C, k,t, e , ei) 

{1 + E\x) + (1 + E(x)[(E k+1 , Y) + hc(k)(E k (x i )E(x j ) + E(x i )E k (x j ))]}, 

so that the contribution of the second term from 06.21) does not exceed 



K 



(C, t, k, m, e , ei ) \\g\\ (1+Ek)9m (l + E k , Y) m - 2 h 3 ^(1 + E(x { ) + E( Xj )){l + E k { Xi ) + E k 

+ (1 + E{ Xi ) + E( Xj ))[{E k+1 , Y) + hc(k)(E k (x i )E(x J ) + E(x l )E k (x J ))]} 2 , 

which again by dividing this sum into two parts, where E(xi) > E(xj) and respectively vice 
versa, reduces to 

k(C, t, fc, m, e , e 1 )||^|| (1+£fc)8m (l + £ fc , Y) m - 2 h 3 ^(1 + + E\ Xi ) 

+ (1 + + E{x j ))[(E'* 1 1 Y) + ^(^(x^^)]} 2 . 

This is again estimated by (16. 3p . It follows now from (16.11) and Proposition 13.11 that 



\T f F - T h F 



t ± lie, 



(l+B'=+3,.)(l+B3,-)(l+-B fe ,-) m_1 hS 



1 (M e h ° s ei (x)) - hn(C,t, k, m, e , ei)||p||( 1+B A,)®m, 



which is the same as (12. 5p . The proof of (12. 6p is quite the same. It only uses Proposition 15.41 
instead of Proposition 15.21 



7 Auxiliary Estimates 

The main technical ingredient in the proof of a weak form of CLT (convergence for fixed times, 
stated in Theorems 12.21 12~~4"|) is given by the following corollary to Theorem 12. 1[ 

Proposition 7.1 Under the assumptions of Theorem \2.1\ let g 2 be a symmetric continuous 
function on X 2 . Then for any k > 1 



sup 

s<t 



E U, 



^(Z fe )-/i s (/i ) 



sup \E(g 2 , (i^ \/io)m = suv\(U%°> s (g 2 , .))(F h ) 



s<t 



s<t 



(7.1) 



does not exceed the expression 



/t(C,t,fc,e ,e 1 )||( ?2 || (1+i;fc) ^ (X 2 ) (l + ( J E; /£ + 3 ,4 + / uo)) 2 I 1 + 



Vh 



1+E k , 
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for any k > 1 under the condition (CI) or (C2) and the expression 



n(C,t, k,eo,ei)\\g 2 \\ c 2 



°(1 + Sfc)^2 (x 2 ) 



1 + {E k+ \ Zl + // )) 3 I 1 + 



Vh 



M 1 



1+E K 



for any k > under the condition ( C3) with a constant k(C, t, k, eo, ei). 
Proof. One has 



E U, 



E to 



to 



Vh 



+ 2E U 



Vh 



Vh 



(7.2) 



The first term can be rewritten as 
1 



h 



E(to (Zt(Zo))® 2 — {ih{Zq)) 



+ m^) {v t {z%) - zt{z%)) + (M^) - ^(4)) ® M4)). 

Under the condition (CI) or (C2) this term can be estimated by 

k{C, r, e , ei) |MI ( i+^(l + (E k+ \ Z*))(l + (£ fc , Z*))(l + (£ 3 ,Z£)) 

<K(C,r,e ,e 1 )||( ?2 || (1+i , fc) « 2 (l + ( J E /c+3 ,4)) 2 , 
due to Theorem 12.11 and (15. 2p . The second term is estimated by 



|to'l(i+E fc )® 2 (l + {E +1 ,Ho + Zq)) 



Zq - Ho 



Vh 



1+E k 



by (12. 3p . and the third term by the obvious combination of these two estimates completing the 
proof for cases (CI) and (C2). The case (C3) is considered analogously. Namely, the first term 
in the representation (17. 2\\ is again estimated by Theorem 12. 11 and to estimate the second term 
one uses ( 15. 15ft instead of ( 12. 3ft and the observation that 



| (to O I < sup 



■t-i 



{i + E k { Xl )y l 



gg2 



xi, x 2 )v{dx 2 ) 



WWm 1 



1+E k 



< 



sup 

Xl,X2 



(1 + E k (x 1 )r\l + E k (x 2 )) 



-i d 2 g 2 



\ v \\m> h 

1 + E k 



< 



\g 2 \\ (c 2,sy 
K (1 + 1 



\\ v M 1 ■ 



dx\dx 2 

Though the estimates of Proposition 17.11 are sufficient to prove Theorem 12.21 in order to 
prove the semigroup convergence from Theorem 12.41 one needs a slightly more general estimate, 
which in turn requires a more general form of LLN, than presented in Theorem 12.11 We shall 
give now these two extensions. 
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Proposition 7.2 The estimates on the r.h.s. of (12. 5p and ( \2.6\i remain valid, if one the l.h.s. 

on takes a more general expression, namely 

sup \T t h (GFH)(Y) - G(Y t )F(Y t )T t h H(Y)\, 

s<t 

where F(Y) is as in Theorem \2. 1\ and both G and H are cylindrical functionals of the form ( 12.121) 



in 



with f G C (R ) and all 4>j, j = 1, ...,n, belonging to Ci +E k(X) and C' Ek (X) respectively 
cases (C1)-(C2) and (C3) (with a constant C depending on the corresponding norms of <pj). 

Proof. As 

T t h (GFH)(Y) - G(Y t )F(Y t )T t h H(Y) 

= T t h (GFH)(Y) - (GFH) (Y t ) + (GF)(Y t )(H(Y t ) - T t h H(Y)), 

it is enough to consider the case without a function H involved. And in this case looking 
through the proof of Theorem 12.11 above one sees that it generalizes straightforwardly to give 
the result required. 

Proposition 7.3 The estimates of Proposition \7.1\ remain valid if instead of (17. ip one takes 
a more general expression 

sup \E[(g 2 , F^(Z%,h ))G(FI?(Z$, // ))] | = sup \[U^ s ((g 2 , .)G)](F Q h )\, (7.3) 

s<t s<t 

where G is as in the previous Proposition (with a constant C again depending on the norms of 
<pj in the representation of G as a cylindrical function of the form (I2.12p ). 

Remark. Let us stress for clarity that U^f ,s ((g2,-)G) in (17. 3p means the result of the 

evolution U^f ,s applied to the function of Y given by (g 2 , Y m )G(Y). 

Proof. It is again obtained by a straightforward generalization of the proof of Proposition 
17.11 given above using Proposition 17.21 instead of Theorem 12.11 

The main technical ingredient in the proof of the functional CLT (stated as Theorems 12.51 
12. 61) is given by the following 



Proposition 7.4 Under condition (C3) for any k > 1/2 

sup E||F s h '(4, /i ) ||? a , < k(C, t, k, e , e x )(l + {E k+ \ Z% + /i )) 3 (l + ||itf ll^o y ). (7.4) 

Proof. The idea is to represent the l.h.s. of ( I7.4p in the form of the l.h.s. of (17. II) with an 
appropriate function g 2 . Using the notation v(x) = f°° v(dy) from the introduction for a finite 
(signed) measure v on R + (and setting v{x) = for x < 0) one has 

= -= / e -**( / u(dx))dy = -= / u(dx) / e-^dy, 

V ZTl Jo Jy V Z7T Jo Jo 

so that for fk(x) = 1 + x k 

T{f k V) = (1 + (i— ) k )Hv) = -7= / v(dx) / (l + y k )e~^dy. 
op V27T Jo Jo 
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Applying (11. lip yields 



\' y W~<^<> o , : - ||/fc^||/r-i(R) — 2^ 



|2 

l(L^(R+))' 



1 

2^ 



— oo 

OO /"CO 



i/(da;) / (1 + y k )e~ ipy dy 



dp 



l+p 2 



o Jo 



6 k {x,y)v(dx)v(dy) 



with 

Clearly 
Moreover 

so that 
and 

Since 



Ok fay) 



CO / PX 



(1 + z k )e~ ipz dz Hi + w k )e ipw dw) -^— r> 
oo \jo Jo J l+p 2 



Ok(x,y)\x=o = O k (x,y)\ y= o = 0. 



dx 



fay) 



\(l + x k )e 



I (1 ■ w k )e ipw dw]-^ 
'o l+P 2 



\-Ji( X ,y)\<c(k)(l + X k )(l+ y k+1 ), 



dxdy 



fay)\ 



d 2 e k 

dx 2 



{l + x k )e' ipx (l + y k )e ipy 

/•CO 

ar,y) = / [kx k ~ l e- ipx 



dp 



l+p 2 
1 + w k )e ipw dw] 



<c(k)(l + x k )(l+y k ) 
dp 



Jo l +P 

r y dn 
[(l + x k )e- ipx / (l + w k )(ip)e ipw dw}—^, 
j Jo l+p 

and using integration by parts in the second term yields also 



dx 2 



fa y)\ < c(k)[(l + x k )(l + y k ) + (1 + ^(l + 



(7.5) 
(7.6) 



Consequently 6 k + 6 k G Cn+^l+i)®2- Therefore, using (17.51) for z/ = F^Zq,^) implies that in 
order to estimate the l.h.s. of (I7.4p one needs to estimate the l.h.s. of (17. ip with g 2 = 9 k given 
by (HU). 

Though a direct application of Proposition 17.11 does not give the result we need, only a 
slight modification is required. Namely, representing (17. ip in form (17. 2p we estimate the first 
term precisely like in the proof of Proposition 17.11 and the second term that now equals 

Vt( z o) -MtOo) ""' 



can be estimated using Proposition 15.81 by 



(T 2 '° V 



7/mII Z MO || 2 



K(c,t, k,eo,ei)(l+E + ,fj, + Z QJU ii,,^-,-- 

Estimating the third term in ( 17.2ft again by the combination of the estimates of the first two 
terms yields ( 17.41) . 
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8 CLT: Proof of Theorems [2J2 - ET6 



Proof of Theorem \2.2\ Recall that we denoted by Uf[ ' r the backward propagator corresponding 
to the process F/ 1 = (Zj* — /j, t )/Vh. By (13.91) . the l.h.s. of (I2.10p can be written as 



sup 

s<t 



(U}i°> s (g,.))(F h )-(U°> s g,F h ) 



sup 

s<t Jo 



[U^(A h T -A T )U^(g,.)]dr(F h ). 



As by (I3~TUD 

(A h T -A T )(U T > s g,-)(Y) 




(U T > s g(y) - U T > s g( Zl ) - U T > s g(z 2 ))K( Zl , z 2 ;dy)Y{dz 1 )Y(dz 2 ) 



Vh 




(U T ' s g(y) - 2U T ' s g(z))K(z, z; dy)(m + VhY)(dz) 



(note that the terms with the second and third variational derivatives in (13.101) vanish here, as 
we apply it to a linear function), the required estimate follows from Proposition 17.11 

Proof of Theorem \2.3l Substituting the function " of form fl 2 . 1 2 [) (with two times 

continuously differentiate ft) with a given initial condition <& r (Y) = '"' n (Y) at t = r in 
the equation F t = —A t F t yields 



9ft df t ,- t . df t 




9 ft 

^ O^riftj' fiy - ~ S z 2 )K(zt, z 2 ; dy)(Y(dzi)nt(dz 2 ) + ^(dzJY ~{dz 2 )) 

j=l 3 




dxjdxi 



4>t, (S y - S Z1 - 5 Z2 )® 2 )K(zi, z 2 ; dy)jit{dz 1 )nt{dz 2 ) 



with ft(xi, x n ) and all its derivatives evaluated at the points Xj = (</>*-, Y) (here and in what 
follows we denote by dot the derivative d/dt with respect to time). This equation is clearly 
satisfied whenever 



f t (xi,...,x n ) = -J2 n (*>$ 
j,k=i 



A) 



d 2 f t 

dxjdxk 



• Xn) 



(8.1) 



and 



{4>){y) - 4>){z) - 4>){w))K{z,w] dy)nt(dw) = -A t $(z) 
with II given by (I2.18p . Consequently 

OU^ r (Y) = $ t (y) = (u^fMu^^Y),...,^^)), 

where W ,r / r = U\[ f r is defined as the resolving operator to the (inverse time) Cauchy problem 
of equation (18.11) (it is well defined as (18. ip is just a spatially invariant second order evolution), 
the resolving operator U t ' r is constructed in Sections 4,5, and 



(8.2) 



n(t, u^tf, u*<ir h ) 
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All statements of Theorem 12.31 follows from the explicit formula (I8.2p . the semigroup property 
of the solution to finite-dimensional equation (18.11) and Propositions I5.4[ 15.71 

Proof of Theorem 2.4 ■ The first statement is obtained by a straightforward modification of 
our proof of Theorem 12.21 above, where one has to use Proposition 17.31 instead of its particular 
case Proposition 17.11 and to note that all terms in formula (13.101) (that unlike the linear case 
now become relevant) depend at most quadratically on Y, because for a function $ of form 

j=i 3 

5 2 $(Y;x,y) = (8.3) 

i,j=l 3 1 

where the derivatives of / are evaluated at the points Xj = (0*-, Y). 

The second statement follows by the usual approximation of a general $ by those given by 
fl2TT2|) with smooth /. 

Proof of Theorem \2.5[ The characteristic function of is 

n 

9L..,tM>->Pn) = Eexp{^(0,,^(4 )/U o))} = U^^.U^^^U^-^^), 

where &j{Y) = exp{ipj(<ftj, Y)}. Let us show that it converges to the characteristic function 
9t u ...,t n (Pu -,Pn) = OU°' t ^ 1 ...OU tn -^-^ n ^OU t --^ n (F ) (8.4) 

of a Gaussian random variable. For n = 1 it follows from Theorem 12.41 For n > 1 one can 
write 

9ti,...,tn(Pl> _ 9t u ...,t n (PU -,Pn) 

n 

= Uff' tl ^ 1 ...U^- 2 ' tj - 1 ^ j -i{U^f j - 1 ' tj - OU t i- 1 > t i)$ j OU t t> t i +1 ...OU tn - ltn <$> n . (8.5) 

3=1 

By Theorem 12.41 we know that for any j = 1, ...,n 

\Uf'~' Pn (Y) = (U h f f ] - 1,t] - OU t i- 1 ' t i)<S> j OU t *> t i+ 1 ...OU tn - ltn § n (Y) 

converge to zero as \fh as h — > uniformly on Y from bounded domains of Ai\ +Ek - We have 
to show that 

•' <I', ,M/,(V) = KU'l>, :*AYr J) (8.6) 

tends to zero, where Ey is of course the expectation with respect to the fluctuation process 
started in Y at time t ^ _ 2 - The last expression can be written as 

E^ldlY^li 2>0 <^ } ^-i^ i )(V 1 )) + E r((l{||y t ,_ 1 || a2 ,o , :))• (8-7) 

( 2>/fc+2 ) ( 2.4+2 ' 

For Y from a bounded subset of (L^ f k 2 )' the second term can be made arbitrary small by 
choosing large enough K due to Proposition 17.41 Due to the natural continuous inclusion 
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C™ k '° C L™f k a , Tn = 1,2 , a > 1/2 one gets by duality a continuous projection (L^)' i— ► 
-Mj fe q C ■M 1 f k _ a for > 1/2, a G (1/2, k). Hence a bounded set in (L\ '° )' is also bounded in 

A^J- fc+2 _ , so that there ^^^(It. J is small of order y/h, implying that the first term in (18.71) 
is small. Consequently expression (I8.6P tends to zero uniformly for Y from bounded domain of 
(L 2 2 '° )', k > 1/2. This implies that all terms in (18.51) tend to zero as h — > 0. 

It remains to check that (18.41) is given by (12 . 1 T[) . which is done by induction in n using 
Theorem 12.31 and an obvious explicit formula 

U l ' r f(x) = expii^pjxj - VjVk \ n(s, <f>$, <j) s k ) ds} 

3=1 j,k=l Jt 

for the solution of the Cauchy problem of the diffusion equation (18.11) with f(x) = exp{i Y^j=i Pj x j}- 
For instance, 

(Oi/*n-i.*»$ n )(y) = (U tn - Utn f n ){U tn - utn <p n ,Y) 

ftn 

= expitpniU^^n, Y)- P l U(S, U'^fa, U^^n) d S } 

Jtn-l 

where f n (x) = exp{ip n x}, and hence 

0£/*»-aA-i($ n _ 2 of7*n-i,*n$ n )(y) 

= exp^p^U^- 2 ^^ + p n U^-»^ n , Y) - pi I " U{s, U s ^(P n , U s > tn (f> n ) ds} 

Jtn-1 

rtn-i 

x exp{- / [pl^Uis, ^""Vn-i, C^^n-i) + 2p n . lPn U(s } U^^, f/ s ^0„)] ds}. 

Jtn-2 



The proof is complete. 
Proof of Theorem \2.6[ 

(i) Notice first that applying Dynkin's formula to the Markov process one finds that for 
a G C 1+E k(X) 

M$(t) = (0, Zt) - (0, Z h ) - f (1^(0, .)){Z h s )ds 

Jo 

is a martingale, since all three terms here are integrable, due to formula (13 .7p and the assumption 
Zq G A4 1+E k+5. Hence (0, F/ 1 ) is a semimartingale and 

M h 

with 



V £(t) = —[^zz) + y o (L,(0,.))(z s ft )rf S -(0,^)] 

is the canonical representation of the semimartingale (0, F^) into the sum of a martingale and 
a predictable process of bounded variation that is also continuous and integrable. (It implies, 
in particular that (0, F/ 1 ) belongs to the class of special semimartingales.) 
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As we know already the convergence of finite dimensional distributions, to prove (i) one has 
to show that the distribution on the Skorohod space of cadlag functions of the semimartingale 
(</>, F/ 1 ) is tight, which according to Aldous-Rebolledo Criterion (see e.g. [8], [30], we cite the 
formulation from [8]) amounts to showing that given a sequence of h n — > as n — > oo and a 
sequence of stopping times r n bounded by a constant T and an arbitrary e > there exist S > 
and n > such that 

sup sup P [\V (n \r n + 9)- V (n \r n )\ > e] < e, 

n>n 0e[O,<5] 



and 



sup sup P [\Q {n \r n + 9)- Q {n \r n )\ > e] < e, 

n>n 6»e[0,5] 



where V n (t) is a shorter notation for V^ n and Q n (t) is the quadratic variation of the martingale 
M^ n (t). Notice that it is enough to show the tightness of (0, F/ 1 ) for a dense subspace of the 
test functions (f>. Thus we can and will consider now only the bounded (f>. 
To get a required estimate for V n (t) observe that by (13. 7p 

jV n (t) = - = ~J J [cj>(y)-2<j>(z)]K(z,z- } dy)Z?(dz) 

+ -^= [ [ ![<j)(y) - <f>{ Zl ) - ( j)(z 2 )]K(z u Z2]dy)[Z^{dz 1 )Z^{dz 2 ) - i^dzMdz*)]. 



2Vh. 

The first term here is clearly uniformly bounded for h — > 0, and the second term can be written 

as 

Iff Jmy)-4 > (z 1 )-4 > (z 2 )]K(z 1 ,z 2 -dy)[Fl l (dz 1 )Z^^ (8.8) 
Applying Doob's maximal inequality to the martingale 

(<^)-(0,F o V f .)){Z h s )ds 
Jo 

in combination with Proposition 17.41 shows that (18. 8p can be made bounded with an arbitrary 
small probability, implying the required estimate for V n (t). 

Let us estimate the quadratic variation by the same arguments as in [16]. Namely, as the 
process (cp, F/ 1 ) for each h is the sum of a differentiable process and a pure jump process, both 
having locally finite variation, its quadratic variation coincides with that of Mr™(t) and is 
known to equal the sum of the squares of the sizes of all its jumps (see e.g. Theorem 26.6 in 
[12]), so that 

Q(n) {t) _ Q(n) (r) = £ ^ ph n _ pKf = * £ (0, Z ^ - Zfc) 2 . 

se[r,t] s<=[r,t] 

As each jump of Z!? is the change of hS x + hS y to h5 x+y one concludes that 

\Q {n \t) - Q {n) (r)\ < hsup \4>\\N t - Nt\ 
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T 



with N t denoting the number of jumps on the interval [0, t]. By the Levy formula for Markov 
chains (see e.g. [3]) the process N t — J* a(Z^) ds is a martingale, where a{Y) denotes the 
intensity of jumps at Y, given by (12. 4p . Hence, using the optional sampling theorem and (12.41) 
implies that 

B(N t -N t ) = eJ a{Z h s ) ds < 3Ch' 1 e (e 1 + e )E(t 
and consequently 

B\Q n (t)-Q n (r)\<3CU\\e (e 1 + e )9 

uniformly for all < t — 9 < t < t. Hence by Chebyshev inequality the required estimate for Q n 
follows. 

(ii) By Theorem 12.51 the limiting process is uniquely defined whenever it exists. Hence 
one only needs to prove the tightness of the family of normalized fluctuations F^. Again due 
to the existence of finite dimensional limits and general convergence theorems (see either a 
result of [26] specially designed to show convergence in Hilbert spaces, or a more general result 
on convergence of a complete separable metric space valued processes in [9] or [8]), to prove 
tightness it is enough to establish the following compact containment condition: for every e > 
and T > there exists K > such that for any h 

P(sup \\F t h \\ {L20 y >K)<e. 
te[o,T] 2J fc 

To this end, let us introduce a regularized square root function R, i.e. R(x) is an infinitely 
smooth increasing function R + \—> R + such that R(x) = y/x for x > 1, and the corresponding 
"regularized norm" functional on (L^ °f )': 

G(Y) = R((Y, Y) (l20 y ) = R((9 k , Y®Y)), 

where 9^ is given by (17. 5p (see Proposition I7.4p . By Dynkin's formula one can conclude that 
the process 



M t = G(F t h ) - G(F h ) - [ l A h a G(Fs) ds 

Jo 



is a martingale whenever all terms in this expression have finite expectations. (Note that we use 
here a more general than usual version of Dynkin's formula with a time dependent generator; 
the reduction of time nonhomogeneous case to the standard situation by including time as 
an additional coordinate of a Markov process under consideration is explained e.g. in [TO]-) 
Expectation of G{F^) is bounded by Proposition 17.41 Moreover, taking into account (18.31) 
and the fact that R^ k \s) = 0{s^- k ) for s > 1, one sees from formulas (I3.10p and (I2.8P that 
A a G(F a ) grows at most quadratically in F a , which again by Proposition 7.4 implies the uniform 
boundedness of the expectation of this term. Applying to M t Doob's maximal inequality yields 
the required compact containment completing the proof of the theorem. 

9 Three lemmas 

We present here three general (not connected to each other) analytic facts used in the main 
body of the paper. Recall that classes C l (M.f(X)) were defined in the introduction. 
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Lemma 9.1 (i) If F G C l (M f (X)) and Y^ G M f (X), then 

F{Y + - F(Y) = I (SHY + s& 0,0 ds. (9.1) 



(ii) If F G C 2 (A / i/ i< /,(X)) or F E C 3 (JAfrf(X)), the following Taylor expansion holds re- 
spectively: 



(a) 



FX + - FQO = (^(1"; 0,O+ /'(l - s)(5 2 F(F + ^; 0, £ ® ds, 

Jo 



(b) fx + - *W = (^(^1 0,O + •, 0, £ ® 

+ IJ\ 1 - sf^Fx + •, •, o, e 3 ) ds. (9.2) 

(Hi) Ift i— > /i 4 G A^/(X) zs continuous in the *-weak topology of Aif(X) and is continuously 
differentiable in the *-weak topology of j\4$(X), then for any F G ^(M.f^X)), (f> <f } 

±F(to) = {6F(fH]-),ik). (9.3) 

Proof, (i) Using the representation 

FX + s(5 x + 8 V )) - F(y) = FX + «4) - FX) + / 5F(F + s5 x + h6 y ; y) dh 

Jo 

for arbitrary points x, y and the uniform continuity of SFX + s $x + hS y ;y) in s, ft, allows to 
deduce from (j!.12p the existence of the limit 

lim -(FX + s(5 x + 5 y )) - FX)) = SFX; x) + 5F(F; y). 

s— »0+ S 

Extending similarly to the arbitrary number of points one obtains (19.11) for £ being an arbitrary 
finite sum of the Dirac measures 8 Xl + ... + 5 Xn . 

Assume now that £ G JAf(X) and — > £ as k — > oo ^-weakly in .M/(X), where all are 
finite sums of Dirac measures. We are going to pass to the limit k — > oo in the equation (19.11) 
written for As F G C(M.f) one has 

F(y + 6) - FX) - F(y + - FX), k-^oo. 

Next, the difference 

I (SF(Y + sC k ;-),Ck)ds- [ (6FX + st;-),0ds 
Jo Jo 

can be written as 

/ (6F(Y + sfa.),Z k -Z)ds+ f (6F{Y + st k ;-)-6FX + st;-),Ods. 
Jo Jo 
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The second term tends to zero, because by our assumption the variational derivation SF maps 
M.f{X) continuously to Cf )00 (X). The first term tends to zero, because — > £ weakly and the 
family of functions SF{Y + s£k] .) is compact in Cf )00 (X) (which is again due to the assumed 
continuity of the derivation SF). 

Statement (ii) is straightforward from the usual Taylor expansion. Turning to (iii) observe 
that 

j t F(fi t ) = hm^(F( f , t+h )-F(fi t )), 
which by (i) and the assumed continuous differentiability can be written as 



lim / ( SF(pt + s(n t+h 
h-+oJ \ 



h 

We want to show that it equals the r.h.s. of (19.31) . We have 



AO; •)> T / fa+rdr ) ds. 



J ^5F(fjLt + s(fjLt +h - m)\ -)>^ J fa+rdr^j ds- 



(SF(^;.),fi t ) 



SF(ii t ;.),^ J fa +T dr - fa j + J \SF(ii t + s{^ t+h - /^); .) - SF(n t ; •), ~ J fa +T dr \ ds. 

The first term here tends to zero as h — > by the weak continuity of fa. The second term tends 
to zero, because the family of measures 

Vt+h ~ l-k = t I fa+rdr 
h Jo 

is bounded and hence compact in the *-weak topology of M^X). 

Lemma 9.2 Suppose S is a compact subset of a linear topological space Y (we are interested 
in the case when Y is a topological dual to a Banach space equipped with its -k-weak topology) 
and Z t is a Markov process on S specified by its Feller semigroup ty t on C(S) with a bounded 
generator A. Let Vt t (z) = (z — £t)/a be a family of linear transformation on Y , where a is a 
positive constant and £ i; t > 0, is a differentiable curve in Y . Let 

®-[0,T\(S) = U te [o,T]^t(5') 

for a T > 0. Then Y t = Q t (Z t ), t G [0,T], is a Markov process in £1\o,t](>3) f or an U T > with 
the dynamics of averages (propagator) 

U'*f(y) = V s J(Y t ) 

given by the formula 

U s > t f(y) = n- 1 ^^ s n t f(y) (9.4) 

for any f G C (Q[o t T\(S)) , t < T, where Q t f(y) = f(&t(y))- Moreover, if such a function f is 
uniformly continuously differentiable in the direction £ t , i.e. if the limit 

lim ~ (f(Q t+T (y)) - f(n t (y))) = -i(y it f)(n t (y)) = V 6 /(fi*(y)) 
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exists and is uniform in f2[o,T](S'), then for all s < t 



jU^f = U s >%f, (9.5) 



where the operator A t is given by the formula 



A t f^n- 1 AQ t f--V i J. (9.6) 

Proof. Formula (I9.4p follows from the definitions of ^ t and Q t - Formulas (19. 5p . (19. 7p follow 
by differentiating (19.41) using the product rule and taking into account that the derivative V/ 
is supposed to be uniform. 

Remark. Similarly, using the identity 

one shows that 

-^-U^f = -KU^f, (9.7) 

holds for s = t. However, to extend this to s < t one needs some additional assumptions on 
the smoothness of the semigroup ty t . 

Lemma 9.3 fl6j / Let Y be a measurable space and the mapping t i— > /j, t from [0,T] to A4(Y) 
is continuously differentiable in the sense of the norm in A4(Y) with a (continuous) derivative 
fit = Vt- Let a t denote a density of fi t with respect to its total variation \fi t \, i.e. the class 
of measurable functions (equivalence is defined as the a.s. equality with respect to the measure 
\/i t \) taking three values —1,0,1 and such that /i 4 = a t \^t\ and = a t \it almost surely with 
respect to \/i t \. Then there exists a measurable function ft{x) on [0, T] x Y such that f t is a 
representative of class a t for any t e [0, T] and 

1 1 A** 1 1 = IImoII + ds f s (y)u s (dy). 
Jo Jy 



We refer for a proof to the Appendix of [16J noting only that f t could be chosen as such a 
representative of a t , which is at the same time a representative of the class of the densities of 
with respect to its total variation measure \u^\, where is a singular part of u t in its Lebesgue 
decomposition with respect to |/^|. 



A On the evolutions with integral generators 

Here we present an analytic study of evolutions with integral generators that are obtained as 
certain perturbations of positivity preserving evolutions. As always, it is assumed that X is a 
locally compact space (though this assumption is used only in Theorem IA.21 other statement 
being valid for arbitrary topological spaces). 
We shall start with the problem 

u t (x) = A t u t (x) = / u t (z)v t (x; dz) — a t (x)u t (x), u r (x) = 4>{x), t>r>0, (A.l) 
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where (j> and at are given measurable functions on X such that at is non-negative and locally 
bounded in t for each x, ut(x, •) is a given family of finite (non-negative) measures on X 
depending measurably on t > 0, x G X, and such that sup te r 0)T i \\vt(%, Oil * s bounded for 
arbitrary T and x. 

Clearly equation (1A.1[) is formally equivalent to the integral equation 

Ut (x) = FAu) t = e - teW ^ rW V(x) + / e-^ t{x) -^ x)) L s u s (x)ds, (A.2) 



where £t(x) = J^a s {x)ds and L t v{x) = J v (z)u t (x, dz). 

We shall look for the solutions of (1A.2I) in the class of functions Ut(x), t > r, that are 
continuous in t (for each x), measurable in x and such that the integral in the expression 
for Lu s is well defined in the Lebesgue sense. Basic obvious observation about (I A. 21) is the 
following: the iterations of the mapping K form flA.2j) are connected with the partial sums 



m 

"(Ct-€r) 



+ ^2 e~^ t -^ ) L Sl ■ ■ ■ L^e-^n-i-^L^e-^i-^dsi ■ ■ ■ ds t 

l = l J r<si<---<si<t 



(where e & designates the operator of multiplication by e &( x )) of the perturbation series 
solution S f ' r = lim^oo to (1A~2D by 

(ijr(0) t = 5K-^+ / e-^^L^-.-L^e-^-^)^^---^. (A.3) 

J r<s m <---<si<t 



Lemma A.l Suppose 

A t iP(x)<aP(x), te[0,T\, (A.4) 
for a strictly positive measurable function if) on X and a constant c = c(T). Then 

(i;) m Wt < U + c(t - r) + ■ ■ ■ + -^c m (t - r) m ^j ifj (A.5) 

for all < r < t < T , and consequently S t ' r ip is well defined as a convergent series for each 
t, x and 

S l ' T iP(x) < e c{t ~ r) iP(x). (A.6) 
Proof. This is given by induction in m. Suppose (jA.50 holds for m. Since (IA.40 implies 

L t ifj{x) < (c + a t {x))if>(x) = (c + £ t (x))if)(x), 

it follows that 

{r^) m+1 {ip)t < e-^ t{x) - &ix)) if){x) 

+ f e-^ tix) -^ x)) {c + i s {x)) (l + c{s-r) + --- + ^-c m {s-r) m ]if){x)ds. 
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Consequently, as 

£ e-^-^i^is - r) 1 ds = i(t - r) L - * J* e -«'-«(s - r) 1 ' 1 ds 
for I > 0, it remains to show that 

m 



\{t - r) 1 - 1 /" r '(-s - ,-)'- V.s 
t. — 1J 



+ E^/ 

z=o Jr 



e -fo-M(s-r) l ds 



< c(t - r) + • • • + 



-c m+1 (t-r^ m+1 



(m + 1)! 

But this holds, because the l.h.s. of this inequality equals 



rn J. r m+l ft 

S TiT(t-r) l + — - / e-&-^(s-r) m ds. 

^ l - m - Jr 



The following corollary plays an important role in the analysis of Section 5. 
Lemma A. 2 Suppose A t ip < cip + <p for positive functions <p and ip and all t e [0, T\. Then 

J r 

Proof. Using (1A.5|) yields 

fyMt < (l + c(t-r))V+ / e-^-^<f)ds, 



{Ilf{^)t < I 1 + c(t - r) + l(t - rf ) V + / e-«^(i + c(a - r))0ds 



+ y e- (5i ~ €s) L s / e~^ s -^(f)dTds, 

etc, and hence 

(/;nv>) t < e c{t ~ r) 



-(&-£. 



+ / e-^~ is) L s f e- {is - &) (f)dTds + --- 

J r J r 



At-r) 



and the proof is completed by noting that 



^ + / ^ ( e ^ ^ + / e'^^^^e - ^"^ ds + ■ ■ ■ 



S*il> = lim < lim (/;) m (V)*- 



The existence of the solutions to (lA.ll) and (|A.2[) can be easily established now. 
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Proposition A.l Under the assumptions of Lemma \A.l\ the following holds. 

(i) For an arbitrary G B^ the perturbation series S t,r <j) = lim m ^oo S^(j) is absolutely 
convergent for allt,x, the function S t,r (f) solves (1A.2I) and represents its minimal solution (i.e. 
S t (f) < u point-wise for any other solution u to (1A.2|) ), and S t,r (f)(x) tends to S T,r (f)(x) as t — > r 
uniformly on any set where both a t and ip are bounded. 

(ii) The family S t,r form a propagator in B^(X) with the norm 

||5'- r |U < c cit ~ r) . (A.7) 

Proof. Applying Lemma I A. II separately to the positive and negative part of <fi one obtains 
the convergence of series S t,r 4> and the estimate (1A.7I) . Clearly S t,r (j) satisfies (1A.2I) and it is 
minimal, as any solution u of this equation satisfies the equation u t = (I^) m (u)t and hence (due 
to (1A.3I) ) also the inequality u t > S, : 



t,r 
m—l 



The continuity of S t,r in t follows from the formula 

+ j (e-&-M - l)e~^-^L s S s ' r <j)ds + e~&-^ L s S s ' r (f) ds 



(A.8) 



for r < t < t. 

At last, once the convergence of the series S t,r is proved, the propagator (or Chapman- 
Kolmogorov) equation (11 .13)) follows from simple standard manipulations with integrals that 
we omit. 

For the application to time non-homogeneous stochastic processes one needs actually equa- 
tion flA.lj) in inverse time, i.e. the problem 



u t (x) = -A t u t (x) = - J u t (z)u t (x; dz) + a t (x)u t (x), u r (x) = <f>(x), < t < r, (A. 9) 
with the corresponding integral equation taking the form 

e^ x) -^ x) L s u s {x)ds. (A.10) 

All the statements of Proposition IA.1I (and their proofs) obviously hold for the perturbation 
series S t,r constructed from flA.101) . with the same estimate (IA.7I) . but with the backward 
propagator equation (11.131) holding for t < s < r with S instead of U. 

To get a strong continuity of S'*' r one usually needs a second bound for A t . In particular, 
the following holds. 

Proposition A. 2 Suppose now that two measurable functions ip2 on X are given both sat- 
isfying flA.4j) and such that (i) < ipi < ip2, (ii) at is bounded on any set where ip2 is bounded, 
(Hi) ipi G B^ 2!00 (X). Then S t,r , t < r (constructed above for (IA.9I) . (IA.10I) ) is a strongly 
continuous family of operators in B^ 2tOQ (X). 

Proof. By Proposition IA.1I S t,r are bounded in B^ 2 (X). Moreover, as S t,r <f) tends to <fi 
uniformly on the sets where ip2 is bounded, it follows that 

for any G B^ 1 (X), and hence also for any G B^ 2!00 (X), since B^ 1 (X) is dense in B^ 200 (X). 
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Theorem A.l Under the assumptions of Proposition I7Q1 assume additionally that ipi-,ip2 are 
continuous, at is a continuous mapping t \— > C^ 2 /f ljOC and L t is a continuous mapping from t 
to bounded operators i— »■ C^ 2)0O . Then is an invariant core for the propagator S t,r in 



the sense that 



5^0-0 
A r m = hm = Iim 



t-*r,t<r r — t s^r,s>r s — T ' 

^Ls^ s <p = 5*'M S 0, ^S s > r <p = -A s S s ' r (j), t<s<r, (A.ll) 
ds ds 

for all G B^ X {X), with all these limit existing in the Banach topology of B^ 2t00 (X) . Moreover, 
Cfa and C^ 2j00 are invariant under S t,r , so that G^ x is an invariant core of the strongly con- 
tinuous propagator S t,r in C^ 2yOQ . In particular, if a t ,L t do not depend on t, then A generates 
a strongly continuous semigroup on C^ 2yOQ with being an invariant core. 

Proof. The differentiability of S t ' r (f)(x) for each x follows from (1A.8I) (better to say its time 
reversal version). Differentiating equation (lA.lOj) one sees directly that S t,r (f) satisfies ( 1A.9j) and 



al required formulas hold point-wise. To show that they hold in the topology of B^ 2t00 one 
needs to show that the operators A t (<f>) are continuous as functions from t to B^ 2t00 for each 
G B^ x . But this follows directly from our continuity assumptions on a t and L t . 

To show that the space is invariant (and this wold obviously imply all other remaining 
statements), we shall approximate S t,r by the evolutions with bounded intensities. Let \ x 
be a measurable function X \— > [0,1] such that Xn( x ) — 1 f° r ^2 0*0 < n and Xn( x ) — 
for ip2{%) > n + 1. Denote u^(x,dz) = x n (a;)i/ t (:r, dz), a™ = Xn a t, and let S^ r (respectively 
A™) denote the propagators constructed as in Proposition IA.2I (respectively the operators from 
(lA.ip ) but with z/™ and a" instead of v t and a t . Then the propagators Sjf converge strongly 
in the Banach space B^ 2t00 to the propagator S ,<,r . One can deduce this fact from a general 
statement on the convergence of propagators (see e.g. [24]), but a direct proof is even simpler. 
Namely, as S t,r and S^ r are uniformly bounded, it is enough to show the convergence for the 
elements of the invariant core B^ x . For such a one has 

(^-5f)(0)= f -fs t > s S% r <f>ds = f S^iAs-A^S^ds, (A.12) 
Jt ds J t 

where ( 1A.11I) was used. As by invariance S^ r (j) G B^, it follows that (A s — A™)S^ r (j) G B^ 2 
and tends to zero in the form of B^ 2 , as n — > oo, and hence the r.h.s. of (IA.12[) tends to zero 
in Bfa, as n — > oo. 

To complete the proof it remains to observe that as the generators of S^f are bounded, the 
corresponding semigroups preserves continuity (as they can be constructed as the convergent 
exponential series). Hence S t)T preserves the continuity as well, as S t,r <p is a (uniform) limit of 
continuous functions. 

Remark. Choosing a t = \\ut(x,-)\\ and -01 = 1 above yield a pure analytic construction of 
a strongly continuous propagator for a non-homogeneous jump type process. A more familiar 
probabilistic approach can be found e.g. in jl] (at least for the homogeneous case). 

For our purposes we need a perturbed equation (IA.9I) . namely the equation 

Ut = -{A t - B t )u t , u r = 0, 0<t<r, (A.13) 
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where B t are bounded operators in C^ 1; and its dual equation on measures, whose weak form 
is 

j t (g,Zt) = ((At-B t )g,£ t ) £ = £, 0<t<r, (A.14) 

i.e. has to hold for some class of test functions g. Motivated by the standard observation that 
formally equation (1A.13j) is equivalent to the integral equation 



u t = S^-J S^ s B s u s d Sl (A.15) 
whose solution u t = U t,r (f) one expects to obtain through the perturbation series 

U^cj) = S^ r <p - [ S t > s B s S s > r ds+ [ S t > Sl B Sl S Sl > S2 B S2 S S2 ' r ds 1 ds 2 + --- , (A.16) 

Jt Jt<si<s 2 <r 

one arrives at the following result. 

Theorem A. 2 Under the assumptions of Theorem ] A. 1\ suppose that ijh{x) —>■ oo as x — > oo 

and that a strongly continuous family of bounded operators B t : C^ 2 i— > is given. Then 
(i) series ( 1A.16I) is absolutely convergent in C^ 2 (X) for any <ft e C^ 2 (X) so that 



lU^Wc^x) < \\S^ r \\c^ ( x)exp{(r-t) sup \\B S \\ C ^ {X) }, 

t<s<r 



and defines a strongly continuous backward propagator U t,r in C^ 2t00 (X) with C^ x being its 
invariant core (so that the analogues of (lA.llj) hold); 

(ii) the operator V r,s = (U s,r )* form a weakly continuous propagator in M.^ 2 yielding a 
unique (weakly continuous) solution to the Cauchy problem (lA.14j) in the sense that it holds for 
allge Cfr{X); 

(Hi) if f is an arbitrary continuous function tending to zero as x —>■ oo, then the operators 
V r,s = (U r ' s )* are strongly continuous in the norm of Ai^p 2 f and solves a strong version of 
f)A.14j) with derivative taken in the norm topology of Ai^f. 

(iv) at last, if a family A", B^ of operators are given satisfying all the above conditions for 
each uj from an interval and such that A^ — B^ depend strongly continuous on uj as opera- 
tors i — ^ C^ 2tOQ , then the corresponding resolving operators U s,r in Cf 2j0C depend strongly 
continuous on uj and their adjoint operators V r,s depend weakly continuous on uj in M.^ 2 . 

Proof, (i) (1A.16j) converges, because B t are bounded. Other statements then follow directly 
from the corresponding facts about S t,r . 

(ii) The operators V r,s are weakly continuous in M.^ 2 (X) just because they are adjoint to 
strongly continuous operators in C^ 2i00 . Next, the analogue of the third equation in flA.lip for 
U t,r is the equation 

±W> r g = U s > r (A r - B r )g 

that holds in C^, 2t00 (X) for any g G according to (i). Passing to the adjoint operators it 
implies 

^(g,V r ' s Y) = ((A r -B r )g,V r > s Y) 
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showing that V r,s yield a solution to (1A.14j) . To show the uniqueness we shall use the method 
for the reduction of the uniqueness problem to the existence of certain solutions of the adjoint 
problem, see e.g. [25] in the Hilbert space setting and time independent generators. Let 
< a < b < r, X[a,b]( s ) be an indicator function of [a, b], and v G C^ 1 {X). As £7 <,r solve 0A.130 
one deduces that the function 



U S ' r X[a,b]( S ) vds 



solves the problem 

d 



-cj) t = -(A t -B t )(j) t + X[a,b]( s ) v , 0r = O, (A.17) 

in the sense that <p t is continuous and satisfies (1A.17I) everywhere with possible exception of two 
points, where its derivative is not continuous. Now, to prove uniqueness for( 1A.14l) it is enough 
to show that its any solution with £o = vanishes. Assume that £t is a weakly continuous 
function in J\4^ 2 (X) such that £o = and ( 1A.14I) holds for all g G C^ x . Integration by parts, 
( 1A.14|) and weak continuity of £ t imply that 



= |[ =0 = f + ((A - B t )<f> t ,Z t )]dt 

Jo 

whenever <fi t has a uniformly bounded derivatives in C^, 2)00 (X) apart from a finite number of 
points. Using (1A.17I) yield now the equation 

J (v,£ t )dt = 0. 

As it holds for arbitrary < a < b < r, v E C^ 1 (X), it implies that ^ = 0. 
(iii) From (1A.14I) it follows that 



(J9,£r)-(j9,£.)= / (( A t ~ B t )g,£ t )dt, < s < r. (A.18) 

J s 

Approximating any g G B$ X (X') by functions from C 1 p l and using the dominated convergence 



one concludes that ( 1A.18I) holds for g G B^ 1 (X). From this one deduces that £ 4 is an absolutely 
continuous function of t in the norm M.^{X). From boundedness of ^ in J^A^ 2 f(X) (that 
follows from weak continuity) and the weak continuity in Jvl 1 p 1 (X) it follows the continuity in 
A4,,p 2 f(X). At last, again from (1A.18I) one concludes that ^ is continuously differentiable in 
M vlf (X). 

(iv) This is straightforward. Namely, one compares U r,s for various to by a formula similar to 
(1A.12I) . This yields the continuous dependence of U r ' s (p on to for G C^ )1 (X). By approximation 
one extends this result to all 6 G Ct?. 
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